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ABSTRACT 


This  thesis  addresses  the  design  and  analysis  of  discrete 
lateral  autopilots  for  application  to  BTT  missiles. 

The  first  part  reviewed  the  classical  design  and  analysis  of 
the  continuous  uncoupled  yaw  and  roll  channels  as  developed  in 
[Ref.  6].   Then,  applying  analog-to-digital  conversion,  the 
corresponding  discrete  autopilots  were  designed  and  analyzed  in 
terms  of  their  transient  responses. 

The  second  part  utilized  modern  control  design  techniques  for 
the  single-input  discrete  lateral  autopilots.   At  first,  assuming 
availability  of  all  states  for  feedback  purposes,  a  discrete 
state-feedback  autopilot  was  obtained.   Next,  since  the  state 
vector  is  not  always  available  to  direct  measurement,  as  estimator 
was  introduced  to  implement  control.   The  state-feedback  and 
estimator  designs  were  analyzed  for  both  lateral  channels  and 
found  to  have  satisfactory  time  responses. 

Finally,  coupling  the  discrete  pitch  and  roll  channel 
autopilots,  a  state-feedback  and  estimator  were  designed  and 
found  to  be  robust. 
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C  .  slope  of  curve  of  pitching  moment  coefficient  (-.;i)  vs 

*  angle-of-attack 
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**Y  yaw  control  incidence  ( JY  ) 


c  , 


<-"y  side  force  coefficient 

CY  slope  of  curve  of  side  force  coefficient  (c"y)  vs 
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p         roll  rate  about   x„ 
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p  roll  acceleration  about  Xg 

POC  Preferred  Orientation  Control 

POPLAR  Pole  Placement  and  Robustness  Design  Program 

q  dynamic  pressure 

q  pitch  rate  about  yB 

r  yaw  angular  rate  about  2g 

rQ  yaw  angular  rate  command  (coordinated  command) 

r  yaw  angular  acceleration  about  2=g 

S  reference  are  for  coefficients 

STT  Skid-to-Turn  (roll  attitude  stabilized) 

u  velocity  component  in  xg  direction 

v  velocity  component  in  y     direction  (assumed  constant) 

V  constant  missile  flight  path  velocity 

V  missile  velocity  vector 

w  velocity  component  in  2g  direction 

W  missile  weight 

x  body-fixed  roll  axis  (along  axis  of  symmetry,  positive 
forward) 

y  body-fixed  pitch  axis  (positive  forward) 

z  body-fixed  yaw  axis  (forms  right-handed  orthogonal 
system  with  x-,  and  y      ) 

1,  achieved  normal  acceleration  in  2   direction 

n  achieved  normal  acceleration  in  y   direction 
w  '8 

\].£  achieved  normal  acceleration  in  2?^  direction 

1 y  achieved  normal  acceleration  in  y       direction 
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'*c  normal  acceleration  command  from  guidance  computer  in 
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degrees  in-f  direction  and  90  degrees  in  y   direction) 


(i  roll  attitude  (zero  degrees  in  -  Zv     direction  and  90 

degrees  in  y   direction) 
'V 

Q  elevation  Euler  Angle  (second  rotation) 

-\l)        azimuth  Euler  angle  (first  rotation  about  Y     ) 

£        pitch  control  incidence  (positive  tail  incidence 
P       produces  negative  pitching  moment) 

So  commanded  pitch  control  incidence 


yaw  control  incidence  (positive  tail  incidence  produces 
'      negative  yawing  moment) 


Cy       commanded  yaw  control  incidence 


g  roll  control  incidence  (positive  tail  incidence 

*  produces  positive  rolling  moment) 

Co  commanded  roll  control  incidence 

SKs  constant  or  equilibrium  angle-of-attack 

•X       angle-of-attack 

angle-of- sideslip 
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I.   INTRODUCTION 

Modern  tactical  missiles  require  increased  stand-off  ranges 
and  need  to  meet  threats  from  highly  maneuverable  air  targets. 
The  high  maneuverability  of  air  targets  has  directed  the  use  of 
defence  missiles  capable  to  develop  higher  lift  accelerations  and 
more  complex  control  laws.   In  order  to  accomplish  the 
requirement  of  large  stand-off  ranges,  propulsion  systems  using 
air-breathing  engines  have  been  studied  and  developed  in  recent 
years.   The  advent  of  air-breathing  engines  has  naturally  led  to 
the  consideration  of  BTT  missiles  in  order  to  minimize  the  inlet 
angle-of- attack. 

The  necessity  of  more  complex  control  laws  has  introduced  the 
application  of  modern  control  and  estimation  theory,  since  more 
complicated  information  of  the  missiles  states  are  needed.   BTT 
controlled  missiles  are  generally  characterized  by  increased 
maneuverability  and  considerable  drag  reduction  over  conventional 
cruciform,  roll  stabilized  STT  controlled  missiles.   Certain 
limitations  in  technology  [Ref.  1]  have  delayed  the  development 
of  BTT  control  systems,  and  consequently  any  progress  in  the  area 
of  BTT  autopilots. 

Major  technological  improvements  during  the  last  decade,  as 
the  availability  of  advanced  digital  computers,  reopened  the 
issue  and  made  BTT  control  feasible  in  spite  the  added  complexity 
of  control  laws  for  the  autopilots.   In  addition,  certain  types 
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of  ramjet  engines  [Ref.  2],  which  are  candidate  propulsion 
systems  for  modern  tactical  mission  requirements  of  range  and 
high  altitude  [Ref. 3],  have  presented  a  need  for  a  missile  control 
technique  to  maintain  effective  inlet  flow.   This  was  the  main 
reason  for  given  further  impetus  to  the  investigation  and 
development  of  BTT  control. 

Despite  the  fact  that  BTT  steering  may  provide  improved 
performance  for  a  missile  system,  there  are  still  unanswered 
questions  concerning  stability  during  homing  phase,  guidance 
performance,  autopilot  guidance  logic  and  subsystem  requirements. 
All  these  questions  have  to  be  investigated  and  properly  answered 
in  order  for  BTT  steering  to  be  considered  as  a  viable  control 
method  for  high  performance  missiles. 

During  the  past  decade  many  missile  programs  [Ref.  3]  were 
initiated  to  improve  the  capability  of  steering   tactical  missile 
via  BTT  control  with  results  that  have  greatly  advanced  the 
understanding  of  the  various  missile  subsystems.   In  the 
autopilot  area  many  different  types  have  been  designed  and 
developed.   All  of  them  force  the  missile  to  roll  or  bank,  so  that 
the  steering  maneuver  occurs  with  the  missile  axis  oriented  in  a 
specific  or  preferred  direction  with  respect  to  the  incoming 
airstream.   This  class  of  autopilots  is  usually  known  as  POC 
autopilots . 

The  main  criterion  for  the  selection  of  a  particular  type  of 
autopilot  is  based  upon  the  guidance,  airframe  and  propulsion 
system  requirements.   Generally,  missiles  with  either  one  or  two 
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planes  of  symmetry  use  a  POC  autopilot  which  forces  the  missile 
to  bank  in  order  to  turn  as  an  aircraft.   If  this  motion  is 
coordinated,  then  the  autopilot  is  referred  to  as  a  CBTT 
autopilot . 

In  the  guidance  area,  radome  aberation  effects  for  frequency 
guidance  are  of  major  concern  [Ref.  31  and  are  being  investigated 
in  great  extent.   Also,  the  interaction  between  BTT  control, 
antenna  stabilization  and  sensor  orientation  are  some  of  the 
additional  concerns  that  have  to  be  properly  addressed.   However, 
simplified  studies  [Ref.  4],  which  neglect  radome  effects  and 
assume  that  the  missile  motion  is  entirely  coordinated,  have 
proven  that  BTT  control  can  provide  acceptable  performance  with 
roll  rates  that  are  not  excessive  for  autopilot  design.   These 
studies  were  made  for  a  medium  range  area  and  long  suppression 
mission,  and  considered  both  high  lift  (i.e.,  planar)  and 
moderate  lift  (i.e.,  cruciform)  airframe  configurations. 

In  order  to  take  full  advantage  of  CBTT  control,  planar 
airframes  have  been  designed  to  increase  the  lifting  capability 
in  -one  direction  without  the  weight  and  drag  penalty  associated 
with  orthogonal  lifting  surfaces  [Ref.  5].   These  airframes  have 
aerodynamic  properties  characterized  by  increased  potential  to 
enhance  CBTT  control. 

The  present  thesis  addresses  the  design  and  analysis  of 
discrete  lateral  autopilots  for  application  to  CBTT  missiles. 

The  first  part  reviewed  the  design  procedure  of  the  two 
individual  lateral  channels  as  developed  by  Arrow  [Ref.  6].   The 
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design  was  performed  using  classical  techniques  and  involved  the 
uncoupled  yaw  and  roll  channels  for  the  elliptical  and  circular 
airframes  respectively.   The  resulted  continuous  open  loop 
designs  were  analyzed  in  terms  of  their  transient  and  frequency 
responses  and  found  to  be  in  accordance  with  the  desired 
requirements  specified  in  [Ref.  6],   Then,  applying  analog-to- 
digital  conversion,  the  corresponding  discrete  lateral  autopilots 
were  obtained  and  analyzed. 

The  second  part  utilized  modern  control  design  methods  to  the 
already  discussed  discrete  single-input  lateral  autopilots.   This 
allowed  comparison  with  the  preceding  classical  design,  and  more 
importantly  established  a  technique  to  extend  some  of  the  results 
to  the  more  general  multivariable  case.   At  first,  assuming 
availability  of  all  states  for  feedback  purposes,  application  of 
the  Ackermann  formula  led  to  a  discrete  state-feedback  designed 
autopilot.   Next,  since  the  state  vector  of  the  state-feedback 
model  is  not  usually  accessible  to  direct  measurement,  an  estimator 
was  introduced  as  an  additional  dynamic  design  in  order  to 
implement  control  to  the  original  system.   The  state-feedback  and 
estimator  designed  autopilots  were  analyzed  for  both  lateral 
channels  and  found  to  have  satisfactory  responses. 

Finally,  coupling  the  discrete  pitch  and  roll  channel 
autopilots,  the  state-feedback  and  estimator  designs  were 
obtained  and  proved  to  be  robust. 

The  analysis  in  all  the  above  cases  was  performed  using  the 
existed  at  Naval  Postgraduate  School  OPTSYS  and  ORACLS  Fortran 
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program  for  the  continuous  and  discrete  systems  respectively. 
Additionally,  the  last  part  that  dealt  with  the  coupled  pitch  and 
roll  channel  autopilot  utilized  the  POPLAR  design  program 
developed  by  Gordon  [Ref.  7]. 
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II.   CLASSICAL  DESIGN  AND  ANALYSIS  OF  LINEAR  UNCOUPLED 

LATERAL  AUTOPILOTS 

A.   GENERAL 

The  initial  phase  in  the  design  of  lateral  CBTT  autopilots 
involved  the  design  and  analysis  of  the  individual  uncoupled 
lateral  channels  (i.e.,  yaw  and  roll)  with  prescribed 
relationships  between  speeds-of-response.   These  relationships 
when  coupled  with  the  corresponding  ones  of  the  longitudinal 
uncoupled  channel  (i.e.,  pitch)  would  meet  the  requirements  of 
the  overall  CBTT  autopilot. 

The  uncoupled  autopilot  design  method  was  classical  and  used 
a  combination  of  frequency  response  and  root  locus  techniques. 
Utilization  of  this  particular  design  method  led  to  the 
achievement  of  practical  bandwidths  (i.e.,  sufficient  high 
frequency  attenuation),  and.  in  turn  provided  the  range  of  required 
missile  body  angular  rates  and  control  motions.   In  addition,  the 
resulting  design  minimized  the  influence  of  aerodynamic 
variations  on  desired  responses.   The  application  of  the 
uncoupled  channels  to  the  whole  CBTT  autopilot  was  accomplished 
by  an  appropriate  choice  of  the  relative  time  constants  of  the 
individual  channels.   In  order  to  achieve  the  desired  maneuver 
plane  acceleration  the  roll  channel  was  designed  to  have  a  time 
constant  of  0.5  seconds.   The  yaw  uncoupled  channel,  which 
follows  the  roll  motion  to  produce  the  required  coordination 
(i.e.,  minimization  of  sideslip  angle),  is  designed  to  have  a 
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more  rapid  response  with  a  time  constant  of  0.39  seconds  for  the 
circular  airframe.   The  detailed  requirements  for  the  classical 
design  of  the  uncoupled  yaw  and  roll  channel  autopilots  are 
presented  in  Appendix  A. 

A  fixed  flight  condition  (i.e.,  constant  Mach  number  and 
altitude)  was  selected  for  this  preliminary  performance  study. 
Fixed  flight  conditions  are  typically  used  in  autopilot  designs 
to  identify  and  cure  critical  areas  of  concern.   When  autopilot 
requirements  are  satisfied  at  fixed  flight  conditions,  then  areas 
of  concern  caused  by  varying  them  are  addressed.   The  selected 
flight  condition  of  60000  feet  altitude  and  Mach  number  3-95 
provided  sufficient  dynamic  pressure,  so  that  the  missile 
maneuvers  resulted  in  large  enough  angles-of-attack  to  exercise 
sideslip  control.   Aerodynamic  data  for  this  particular  flight 
condition  are  provided  in  Appendix  A. 

The  aerodynamic  models  developed  for  stability  studies  in  the 
frequency  domain  were  linearized  about  a  trim  angle-of-attack  for 
both  lateral  channels.   The  following  three  assumptions  were 
made : 

1.  The  plane  Xq     -     ■£„    of  Figure  2.1  was  the  maneuver  plane. 

2.  The  missile  was  trimmed  in  pitch  (i.e.,  My=  0,  at  fixed 
values  of  X  ,  q,  and  c_. 

Rather  than  use  the  assumption  that  the  missile  roll  rate  is 
approximately  zero  as  it  is  normally  done  for  the  roll  stabilized 
STT  control,  the  following  assumption  was  made  for  BTT: 
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3.   The  missile  roll  rate  was  constant. 
Linearized  aerodynamic  derivatives  are  given  in  Appendix  C. 

In  this  chapter  the  analysis  of  both  continuous  and  discrete 
uncoupled  lateral  channels  of  a  BTT  autopilot  was  based  on  the 
transient  and  frequency  responses  of  maneuver  plane  accelerations, 
body  angular  rates  and  tail  incidence  angles.   A  general  block 
diagram  of  a  BTT  autopilot  with  all  its  channels  is  shown  in 
Figure  2.2.   Inertial  acceleration  command  were  applied  in  polar 
coordinates  (i.e.,  magnitude  of  the  command  \\      applied  to  pitch 
and  the  direction  $.     to  roll  autopilot).   The  yaw  autopilot  was 
slaved  to  the  roll  autopilot  in  order  to  minimize  the  yaw  and 
roll  motions.   Achieved  maneuver  plane  accelerations  in 
rectangular  coordinates  (i.e.,  r)  *  an^  i]v  )  were  determined  by 
resolving  achieved  body-fixed  accelerations  (i.e.,  J]  ,  and  y\y    ) 
through  missile  roll  rate     (i.e.,  Euler  angles  0  and  U)  were 
assumed  to  be  sufficiently  small). 

B.   AIRFRAME  CONFIGURATIONS 

The  two  airframe  configurations  studied  in  this  work  were 
taken  from  [Ref.  6]  and  are  shown  in  Figure  2.3  and  Figure  2.4. 
Although  the  configuration  in  Figure  2.3  reveals  a  body  of 
circular  cross  section  and  that  of  Figure  2.4  an  elliptical  one, 
both  airframes  have  the  same  cross  sectional  area  distribution. 
In  specific,  the  circular  cross  sectional  body  has  a  closure 
ratio  Abase/Amax  of  0.69  with  Amax  occurring  at  68%  missile  body, 
whereas  the  elliptical  airframe  has  a  3:1  cross  section. 
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Both  airframe  configurations  are  tail-controlled  using  four 
identical  control  surfaces  which  are  located  flush  with  the  body 
base  with  a  +30°  dihedral.   In  the  case  of  the  elliptical  body, 
the  hinge  line  was  skewed  such  that  a  10°  control  deflection 
measured  at  the  body-tail  juncture  had  a  resultant  7.04°  surface 
deflection.   Thus,  the  aerodynamic  control  effectiveness  in  terms 
of  deflection  measured  at  the  body-tail  is  lower  for  the 
elliptical  airframe  although  it  is  nearly  the  same  in  terms  of 
resultant  surface  deflection. 

The  total  span  of  the  mono-wings  is  the  same  for  each 
configuration,  which  results  in  larger  wing  area  for  the  circular 
airframe.   The  wing  area  and  span  for  the  circular  airframe  were 
chosen  as  typical  of  current  maneuvering  missiles.   The  wing  for 
the  elliptical  concept  was  determined  by  projecting  the 
elliptical  body  on  the  circular  body-wing  planform.   The  resultant 
when  exposed  wing  planform  became  the  wing  for  the  elliptical 
body . 

Comparison  of  the  elliptical  airframe  with  the  corresponding 
circular  indicates  the  following: 

1.  About  30%  more  normal  force  that  is  nearly  independent  of 
angle-of-attack  can  be  achieved  at  supersonic  speeds. 

2.  Values  of  longitudinal  stability  parameter  Cm  are  more 
positive,  and  with  more  pronounced  nonlinearities  in 
pitching  moment  at  subsonic  speeds. 

3.  Levels  of  directional  stability  are  increased  and  more 
compatible  with  levels  of  longitudinal  stability. 
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4.   More  yaw  control  is  available  although  suitable  locations 
for  tails  on  the  body  are  more  limited  because  of  the 
geometry  of  the  elliptical  airframe. 
The  two  airframe  configurations  were  sized  to  provide 
realistic  geometric  and  mass  properties.   The  details  are 
presented  in  Appendix  D. 

C.   UNCOUPLED  YAW  CHANNEL  AUTOPILOT  FOR  ELLIPTICAL  AIRFRAME 

The  purpose  of  the  uncoupled  yaw  channel  autopilot  of  a  CBTT 

missile  is  to  minimize  the  sideslip  angle  (B),  or  provide 

coordinate  motion  between  the  yaw  and  roll  channels.   The  easiest 

way  to  accomplish  this  is  by  designing  the  uncoupled  yaw  channel 

(i.e.,  roll  and  pitch  dynamic  effects  neglected)  as  a  regulator 

(i.e.,  no  guidance  command  and  with  rate  and  acceleration 

feedback)  to  help  minimize  the  sideslip  angle. 

A  block  diagram  of  the  uncoupled  yaw  channel  is  shown  in 

Figure  2.5.   In  this  diagram  both  the  aerodynamic  model  and  yaw 

control  law  are  involved.   The  normal  acceleration  ( })v  )  is  not 

used  to  command  the  CBTT  autopilot.   Instead,  it  is  used  for  the 

design  and  analysis  of  the  uncoupled  channel.   The  command  used  by 

the  coupled  system  is  shown  in  dashed  lines  and  is  a  yaw  angular 

rate  command  (rQ).   The  yaw  control  law  shown  in  Figure  2.6  [Ref. 

6]  is  governed  by  missile  body  angular  rate  (r)  and  yaw  normal 

acceleration  (Hy).   At  the  flight  condition  of  interest  (i.e., 

60  kft  altitude  and  Mach  number  3.95)  the  yaw  control  law 

determines  the  required  command  (  <rv  )  to  an  actuator  which  is 
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approximated  as  a  first-order  lag  at  30  Hz.   The  rate  compensator 
computes  the  high  frequency  attenuation  and  is  used  to  minimize 
aerodynamic  variations  on  the  quality  of  the  regulator.   On  the 
other  hand,  the  acceleration  compensator  measures  the  acceleration 
bandwidth  via  the  time  constant  of  the  acceleration  response  of  }) 
1 .   Transfer  Functions  of  Aerodynamic  Model 

The  aerodynamic  transfer  functions  of  the  uncoupled  yaw 
channel  autopilot  are: 


i/ 


z  _  K(-AitU)       tc-Ai  +  U)    5fJ 


r 


C 


(deg/sec/deg) 


(II. C. 1-1) 


Y) 


r 


-4bf 6c 


c 


6 


4Ef  6c 


2 

S  V  1 


—  (g's/deg) 


(II. C. 1-2) 


where : 


vv 


r. 


(II. C. 1-3) 


v^ 


r 


(II. C. 1-4) 


-      6"7  5J  4  £  c\ 

t  h 


(II. C. 1-5) 
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c  = 


K* 


(57.5)  q,$A 


I 


zz 


/B54- 
V 


(II. C. 1-6) 


(II. C. 1-7) 


Substituting  the  values  of  aerodynamic  data  (Table  VII), 

linearized  aerodynamic  derivatives  (Table  VIII),  and  geometric 

and  mass  properties  (Table  IX),  equations  (II. C.  1-3)  through 
(II .C. 1-7)  become: 

A  =  -0. 1351  (II. C 1-8) 

B  =  0.0436  (II. C. 1-9) 

C  =  -17.2748  (II. C. 1-10) 

i  =  34.5495  (II. C. 1-1 1 ) 

K  =  0.4823  (II. C. 1-12) 


Introducing  the  above  equations  (II. C. 1-8)  through 
(II. C.  1-12)  to  (II. C. 1-1)  and  (II. C.  1-2),  the  aerodynamic 
transfer  functions  of  the  uncoupled  yaw  channel  for  the  circular 
airframe  and  zero  angle-of-attack  can  be  obtained. 

a.   Transfer  Function  of  Yaw  Angular  Rate: 

'■So"**') 

(4.9     / 

(II. C. 1-13) 


Z 


o- 


r 


CU5 


0.  0748 


Js5n 


I 


■4^S 


5  fi 


h 


b.   Transfer  Function  of  Yaw  Normal  Acceleration: 


7J6    (K4isH)  (tt*") 


4,5 


5f1 


■4.25 


5  f  i 


(II. C. 1-14) 
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Rearranging  the  above  transfer  functions  in  terms  of  the 
variable  pairs  (r,  Gy )    and  (n  ,  or  )  respectively,  they  become: 

rs% 0. 04^4 ZS-1S.Z2I 5  r.--  56. €6 c^s-X, 74-22 ^T  ( II. C.  1-15) 

nYsfao4^4aYs-ig.^8ifnT5?.4^^sSTi'iaooisoys- i^4.6soi  5T     cn.c.i-i6) 


Applying  inverse  Laplace  transformation,  the  following 
set  of  linear  differential  equations  can  be  obtained: 

Cf0.04?4-£  -18.Z8i$c=- 56.66^   -^  74^2  lr  (II  c  1_1?) 

nYfO.C4f4nT-1g,^SJ5fir  =  24^5^r-O.OOIs5Y-li4.  6S01  £Y   (n.c.1-18) 

Both  equations  (II. C.  1-17)  and  (II.C.1-18)  form  a  second- 
order  system  of  linear  differential  equations  in  which  the 
forcing  function  involves  derivative  terms.   Using  rules  of 
state-space  representation  [Ref.  8]  the  following  equations  are 
obtained : 

X,  =  y^  -  36.66  £Y  (II. C. 1-19) 

^  =  1S.^15X-  0.04*4*2-  0,  95>\2Zr  ( II.  C.  1-20) 

*l  r  %  -O.U54  ^y  (II. C. 1-21) 

i>=  19.^15^,-0.0^4  ^-^SO.^SSS^y  (II. C.  1-22) 

2  =  x'i  ( II.  C.  1-23) 

z  '    XJ.~  ^6>66  c^  (II. C. 1-24) 

>1T  =  *,  +1.46+S  >Xy  ( II.  C.  1-25) 
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2 .   Equations  of  Yaw  Control  Law  and  Actuator 

a.  Acceleration  Compensator  Equation 

The   acceleration   compensator   equation   obtained    from 

Figure    2.6    is: 

0.3194-6      r w  s  fTT   r  p   i ) 

Rearranging  and  applying  inverse  Laplace  transforma- 
tion, (II. C. 2-1)  turns  into  the  following  linear  differential 
equation : 
/• -5/ +1.5^3  fly -1.5973  HYc  (ii.  c.  2-2) 

Substituting  equation  (II. C.  1-25)  into  (II. C. 2-2), 
the  last  becomes: 
y  ■=  1, 5975 Z^  -  Sy\  5.  ^354- Zy  -  }.S?73  n^  (II .  C .  2-3  ) 

b.  Rate  Compensator  Equation 

The  rate  compensator  equation  also  obtained  from 
Figure  2.6  is: 

4.35   f4rS+i)  (II. C. 2-4) 

oy     =  — (ytt) 

.5 

Rearranging  and  applying  inverse  Laplace  transforma- 
tion to  (II. C. 2-5),  it  turns  into: 
cy    =  Q,4-%Sy  +  QA8S  i  f  4.S5y+-  £.S5t  (II. C. 2-5) 

Substituting  equations  (II. C.  1-23),  (II. C. 1-24)  and 
(II. C. 2-3)  into  the  above,  it  becomes: 
*r^=  4.  Six,  i-C.+eSXJ  +  o;iH7zil2AiSy-'\5.<*U4-er'O.Tr+7  JlY    (li.c.2-6) 
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Actuator  Equation 

The  actuator  equation  obtained  from  Figure  2.5  is 


1 


Y 


1 


i8&A 


S  +  l 


r. 


(II. C. 2-7) 


Rearranging  and  applying  inverse  Laplace  transforma- 
tion, the  last  equation  turns  into: 
3Yr'ie8.^cy  -!BB,4or  (II. C. 2-8) 

3 .   Design  Approach  and  Analysis  of  Continuous  System 
Utilizing  state-space  representation,  the  equations 
(II. C. 1-19)  through  (II. C. 1-22),  (II. C. 2-2)  (II. C. 2-6)  and 
(II. C. 2-8)  can  be  modeled  in  a  seventh-order  system  of  the  form 
x=Fx+Gu.   The  continuous  plant  system  and  input  matrixes  F  and  G 
are  shown  in  Table  I,  and  the  state  vector  is: 
K 

zi 

y 
\ 

I 

where  the  state  variables  are: 


(II. C. 3-1  ) 


x1  >  x2 
z  -|  ,  Z£ 


y 


I 


yaw  angular  rates 

yaw  normal  accelerations 

output  of  acceleration  compensator 

input  command  in  the  actuator 

yaw  tail  incidence 
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Executing  the  OPTSYS  program,  using  an  input  step 
function  which  represents  "1  gee  command"  at  zero  trim  angle-of- 
attack,  the  pole-zero  and  time  and  frequency  response  plots  are 
obtained . 

The  pole-zero  plot  of  Figure  2.7  indicates  that  the 
continuous  open  loop  system  is  stable,  since  the  s-plane  poles 


are : 


S^-174.376 


5^-6.104-^  +j  IC.6S96 
S5  =  -  6AQ4X  -j  10.6596 

S4  =  -0>02/747s +j  4;09?l<? 
Ss  r-O.C2l7475-j'4C59l^ 

5 7  -  -X .  92.  96  -  j  -?  c/^$ 


(II. C. 3-2) 
(II. C. 3-3) 
(II. C. 3-4) 

(II. C. 3-5) 
(II. C. 3-6) 
(II. C. 3-7) 
(II. C. 3-8) 


The  time  response  plots  of  the  yaw  normal  acceleration,  angular 
rate  and  tail  incidence  are  shown  in  Figures  2.8  through  2.10. 
In  particular  the  yaw  normal  acceleration  time  response  plot  has 
a  0.39  seconds  time  constant,  7%   overshoot  and  a  steady-state 
error  of  0.018.   These  results  are  in  accordance  with  the 
requirements  referred  in  Appendix  A,  that  is  a  time  constant  of 
0.4  seconds,  overshoot  less  than  10%  and  steady-state  error  not 
necessarily  equal  to  zero.   All  the  above  three  time  response 
plots  are  identical  with  those  presented  in  [Ref.  6]. 

Figures  2.11  through  2.16  show  the  frequency  response 
plots  of  the  yaw  normal  acceleration,  angular  rate  and  tail 
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incidence,  from  which  the  phase  crossover  frequencies  and  gain 
margins  of  Table  II  can  be  obtained.   The  positive  gain  and  phase 
margins  of  the  open  loop  system  ensure  the  relative  stability  of 
the  closed  loop  (controlled)  system. 

4 .   Design  Approach  and  Analysis  of  Discrete  System 

Utilizing  analog-to-digital  conversion  by  the  aid  of 
ORACLS  program  and  for  a  sample  period  of  0.0125  seconds,  a 
seventh-order  discrete  system  of  the  form  x(k+1  )=Ax(k)+Bu(k)  is 
obtained.   The  discrete  pant  system  and  input  matrices  A  and  B 
are  shown  in  Table  III. 

The  pole-zero  plot  of  Figure  2.17  indicates  that  the 
discrete  open  loop  system  is  also  stable,  since  the  z-plane 
poles  are: 

Z{-  O.I  1 5076  (II.C.4-1 ) 

^=C.<?!&554-+jai;U8<a  (II. C. 4-2) 

Z3--CWI&3S4-J  a  1^86-2     .  (II.C.4-3) 

Z4  ,  C.  9984U  f  J  O.O*\0r88  ( II .  C .  4-4  ) 

zs  --  M98+U  -J  a0SlO7&8  ( II .  C  .  4-5  ) 

*+0,9S$S6S*jaoS6IS46  (II.C.4-6) 

Z7--a^556^0.O*n346  UI.C.4-7) 

The  time  response  plots  of  the  yaw  normal  acceleration, 
angular  rate  and  tail  incidence  for  the  discrete  uncoupled  yaw 
channel  are  presented  in  Figures  2.18  through  2.20.   A  close 
observation  of  these  plots  indicates  that  they  are  identical  with 
those  of  the  continuous  classical  system  found  in  the  previous 
section  . 
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Figure    2.7 


Pole-Zero  Plot;  Uncoupled  Yaw  Channel  Autopilot; 
Classical  Design  Continuous  Open  Loop  System; 
Elliptical  Airframe 
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Figure    2.8 


Yaw  Normal  Acceleration  vs  Time;  Uncoupled  Yaw 
Channel  Autopilot;  Classical  Design;  Continuous 
Open  Loop  System;  Elliptical  Airframe 
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Figure  2.9 


Yaw  Angular  Rate  vs  Time;  Uncoupled  Yaw  Channel 
Autopilot;  Classical  Design;  Continuous  Open 
Loop  System;  Elliptical  Airframe 
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Figure    2.10 


Yaw  Tail  Incidence  vs  Time;  Uncoupled  Yaw  Channel 
Autopilot;  Classical  Design;  Continuous  Open 
Loop  System;  Elliptical  Airframe 
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Figure  2. 1  1 


Yaw  Normal  Acceleration-Gain  vs  Frequency;  Uncoupled 
Yaw  Channel  Autopilot;  Classical  Design;  Continuous 
Open  Loop  System;  Elliptical  Airframe 
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Figure  2.12 


Yaw  Normal  Acceleration-Phase  vs  Frequency; 
Uncoupled  Yaw  Channel  Autopilot;  Classical 
Design;  Continuous  Open  Loop  System;  Elliptical 
Airframe 
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Figure  2.13 


Yaw  Angular  Rate-Gain  vs  Frequency;  Uncoupled 
Yaw  Channel  Autopilot;  Classical  Design; 
Continuous  Open  Loop  System;  Elliptical 
Airframe 
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Figure  2.14 


Yaw  Angular  Rate-Phase  vs  Frequency;  Uncoupled 
Yaw  Channel  Autopilot;  Classical  Design; 
Continuous  Open  Loop  System;  Elliptical 
Airframe 
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Figure  2.15 


Yaw  Tail  Incidence-Gain  vs  Frequency;  Uncoupled 
Yaw  Channel  Autopilot;  Classical  Design; 
Continuous  Open  Loop  System;  Elliptical 
Airframe 
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Figure  2.16 


Yaw  Tail  Incidence-Phase  vs  Frequency;  Uncoupled 
Yaw  Channel  Autopilot;  Classical  Design; 
Continuous  Open  Loop  System;  Elliptical 
Airframe 
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D.   UNCOUPLED  ROLL  CHANNEL  AUTOPILOT  FOR  CIRCULAR  AIRFRAME 

The  uncoupled  roll  channel  autopilot  of  a  CBTT  missile  is 
commanded  to  roll  the  missile  so  as  to  put  the  preferred  maneuver 
direction  in  the  direction  of  the  guidance  acceleration  command. 
The  desired  maneuver  plane  acceleration  should  be  attained  as 
rapidly  as  the  achieved  body-fixed  pitch  acceleration.   To 
accomplish  this,  the  uncoupled  roll  channel  autopilot  (i.e.,  yaw 
and  roll  dynamic  effects  neglected)  was  designed  to  have  the  roll 
angle  time  constant  equal  to  the  time  constant  of  the  normal 
acceleration  achieved  by  the  uncoupled  pitch  channel  autopilot. 

A  block  diagram  of  the  uncoupled  roll  channel  is  shown  in 
Figure  2.21.   In  this  diagram  both  the  aerodynamic  model  and  roll 
control  law  are  involved.   The  roll  control  law  shown  in  Figure 
2.22  is  commanded  by  roll  angle  (tyQ)    and  governed  by  roll  angular 
rate  (p)  and  roll  angle  ((f)). 

The  design  and  analysis  of  the  uncoupled  roll  channel 
autopilot  was  performed  in  this  section  for  the  stable  at  zero 
angle-of-attack  circular  airframe. 

1  .   Transfer  Functions  of  Aerodynamic  Model 

The  aerodynamic  transfer  functions  of  the  uncoupled  roll 
channel  autopilot  obtained  from  Figure  2.21  are: 

a.   Transfer  Function  of  Roll  Angular  Rate: 

i    4  $d  ctto  c 


t  ~^T      T         "fc  (II. D. 1-1) 
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Figure  2.17 


Pole-Zero  Plot;  Uncoupled  Yaw  Channel  Autopilot; 
Classical  Design;  Discrete  Open  Loop  System; 
Elliptical  Airframe 
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Yaw  Normal  Acceleration  vs  Time;  Uncoupled  Yaw 
Channel  Autopilot;  Classical  Design;  Discrete 
Open  Loop  System;  Elliptical  Airframe 
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Figure    2.19 


Yaw  Angular  Rate  vs  Time;  Uncoupled  Yaw 
Channel  Autopilot;  Classical  Design; 
Discrete  Open  Loop  System;  Elliptical 
Airframe 
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Figure    2.20 


Yaw  Tail  Incidence  vs  Time;  Uncoupled  Yaw 
Channel  Autopilot;  Classical  Design; 
Discrete  Open  Loop  System;  Elliptical 
Airframe 
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b.   Transfer  Function  of  Roll  Angle: 
^   =  —   p  (II. D. 1-2) 

Substituting  the  values  of  aerodynamic  data  (Table  VII), 
linearized  aerodynamic  derivatives  (Table  VIII)  and  geometric  and 
mass  properties  (Table  IX),  equation  (II. D. 1-1)  becomes 

ps  =  8.03^62  c^  (II. D.  1-3) 

Applying  inverse  Laplace  transformation  to  equations 
(II. D. 1-2)  and  (II.  D.  1-3),  the  following  set  of  linear 
differential  equations  is  obtained: 

p  =  8.03462  o,^  (II. D.  1-4) 

<P  =  p  (II. D.  1-5) 

2 .   Equations  of  Roll  Control  Law  and  Actuator 

a.  Roll  Angle  Compensator  Equation 

The  roll  angle  compensator  equation  obtained  from 
Figure  2.22  is: 

2  Z 

X  =   j-  ;  C<Pc-<!>)  (II. D. 2-1) 

Rearranging  and  applying  inverse  Laplace 
transformation  (II. D. 2-1)  turns  into  the  following  linear 
differential  equation: 

X  =  -17.6$-  8X  +  17.6$c  (II.  D.  2-2) 

b.  Rate  Compensator  Equation 

The  rate  compensator  equation  also  obtained  from 
Figure  2.22  is: 

f      (-f-fij)       ^X  V  (II. D. 2-3) 

60 


o 

c 
ca 


o 


CD 

d> 

"O 

cr 
to 


o 

I/) 
CO 


T3 
to 

u 

<o 


CO 


^3 
O 


O 
U 


CM 

cvi 

CD 

!-. 

3 

■H 


g  §  I 

K  o" 
u 


■o 

03 


CJ 

-e. 


61 


■o 

a: 


u 

cc 


t 

. 

«-" 

+ 

*- 

o 

+ 

l/i 

p 

«lr- 

^ 

o 


„ , 

' — 

+ 

1  ID 

t/> 

- 

— ■■ 

f— 

LO 

CO 

+ 

o 

irt  |  lo 

IT) 

^~ 

o 

"*"*' 

o 

T3 


CN 

+ 

CM 

">|00 

-V    1 

^ 

03 

CO 

o 

3 

*-» 

CJ 

L. 

— 

CJ 

UJ 

•^r 

CN 

r^ 

o 

»— 

i- 

•^r 

II 

ii 

X 

^ 

o 


T3 


2 

r3 


O 
5- 

c 
o 


o 

cc 


(VI 
CNJ 

o 

3 
iO 

•  H 


u  "2 


62 


Rearranging  and  applying  inverse  Laplace  transforma- 
tion to  (II. D.  2-3),  it  turns  into: 
Y  +  5Y  =  0.0502X  -  0.0502p  +  0.75325X  -  0.75325p       (II. D. 2-4) 

Substituting  equations  (II. D. 1-4)  and  (II. D. 2-2)  into 
(II. D. 2-4),  the  last  becomes: 

Y=-0.  75  32p-0.88352<4>+0.  35  165X-5Y-0.  403  34  ^+0.88352%    (II.D.2-5) 
c.   Pseudo-Differential  Equation 

The  Dseudo-dif f erential  eauation  obtained  from  Fieurfi 


2.22  is: 

0.  C/30^5  S 

6 


(II. D. 2-6) 


Rearranging  and  applying  inverse  Laplace 
transformation  to  (II. D. 2-6),  it  turns  into: 
X1  =  -6X-,  +  0.078l98p  (II. D. 2-7) 

Substituting  equation  (II. D. 1-4)  into  (II.  D.  2-7),   the 
last  becomes: 

X1  =  -  6X1  +  0.628291  <5"R  (II. D. 2-8) 

d.   Equation  of  Actuator  Compensator 

The  equation  of  actuator  compensator  obtained  from 
Figure  2.22  is: 

s  no    r  Y  v  \ 

°K~  ~i C    l-XJ  (II. D. 2-9) 

Rearranging  and  applying  inverse  Laplace  transforma- 
tion to  (II. D. 2-9),  it  turns  into: 
SR   =-15^£  +0.  13636Y-0.  13636X1  +  15Y-15X1  (II.  D.  2-10) 
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Substituting  equations  (II.  D.  2-5)  and  (II.  D.  2-8)  into 
(II. D. 2-10),  the  last  becomes: 

(III. D. 2-11 ) 
e.   Actuator  Equation 

The  actuator  equation  obtained  from  Figure  2.21  is: 


*   = 


£715 


££A 


5  4-1 


(II. D. 2-12) 


Rearranging  and  applying  inverse  Laplace 
transformation,  the  last  equation  turns  into: 


^  =  10795.32^  -188.4  o^ 


(II. D. 2-13) 


3 .   Design  Approach  and  Analysis  of  Continuous  System 
Utilizing  state-space  representation,  the  equations 
(II. D. 1-4),  (II. D. 1-5),  (II. D. 2-2),  (II. D. 2-5),  (II. D. 2-8), 
(II. D. 2-11),  and  (II. D. 2-13)  can  be  modeled  in  a  seventh-order 
system  of  the  form  x  =  Fx+G-u.   The  continuous  plant  system  and 
input  matrices  F  and  G  are  shown  in  Table  IV,  and  the  state 
vector  is: 


Xc'K]  = 


4> 

X 

Y 

X. 

L  *«  I 


(II. D. 3-1 ) 
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where  the  state   variables  are 


P 

♦ 

X 

Y 
X 
C 


K 


roll  angular  rate 

roll  angle 

output  of  roll  angle  compensator 

output  or  rate  compensator 

output  of  pseudo-differentiator 

input  command  in  the  actuator 

roll  tail  incidence 

Executing  the  OPTSYS  program,  using  an  input  step 
function  which  represents  "1  gee  command"  at  zero  trim  angle-of- 
attack,  the  pole-zero  and  time  and  frequency  response  plots  are 
obtained . 

The  pole-zero  plot  of  Figure  2.23  indicates  that  the 
continuous  open  loop  system  is  stable,  since  the  s-plane  poles 


are 

5 
i 


V 


=  -174.785   (roll  angular  rate) 
=  -9.25097+J28.4098  (roll  angle) 


(II. D. 3-2) 
(II. D. 3-3) 


<.' 


5 


Ss 


c 


S 


and 


=  -9.25097-j28.4098  (output  of  roll  angle  compensator) 

(II. D. 3-4) 

=  -2 . 46608+ j2 . 71 1 52  (output  of  rate  compensator  network) 

(II. D. 3-5) 

=  -2.46608-j2.71 152  (output  of  pseudo-differentiator) 

(II. D. 3-6) 

=  -8.98209  (input  command  in  the  actuator)         (II. D. 3-7) 

=  -5.19932  (roll  tail  incidence)  (II. D. 3-8) 

The  time  response  plots  of  the  roll  angle,  angular  rate 

tail  incidence  are  shown  in  Figure  2.24  through  2.26.   In 
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particular,  the  roll  angle  time  response  plot  has  a  0.55  seconds 
time  constant,  3%    overshoot  and  a  steady-state  equal  to  zero. 
These  results  are  in  accordance  with  the  requirements  referred  in 
Appendix  A,  that  is  a  time  constant  of  0.5  seconds,  overshoot 
less  than  10%  and  zero  steady-state  roll  angle  error.   All  the 
above  three  time  response  plots  are  identical  with  those 
presented  in  [Ref.  6], 

Figures  2.27  through  2.32  show  the  frequency  response 
plots  of  the  roll  angle,  angular  rate  and  tail  incidence,  from 
which  the  phase  crossover  frequencies  and  gain  margins  of  Table 
V  can  be  obtained.   The  positive  gain  and  phase  margins  of  the 
open  loop  system  ensure  the  relative  stability  of  the  closed  loop 
(controlled)  system. 
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Figure  2.23 


Pole-Zero  Plot;  Uncoupled  Roll  Channel  Autopilot; 
Classical  Design;  Continuous  Open  Loop  System; 
Circular  Airframe 
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Figure    2.2M 


Roll  Angle  vs  Time;  Uncoupled  Roll  Channel 
Autopilot;  Classical  Design;  Continuous  Open 
Loop  System;  Circular  Airframe 
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Figure    2.25 


Roll  Angular  Rate  vs  Time;  Uncoupled  Roll  Channel 
Autopilot;  Classical  Design;  Continuous  Open  Loop 
System;  Circular  Airframe 
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Figure    2.26 


Roll  Tail  Incidence  vs  Time; 
Autopilot;  Classical  Design; 
System;  Circular  Airframe 


Uncoupled  Roll  Channel 
Continuous  Open  Loop 
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Figure  2.27 


Roll  Angle-Gain  vs  Frequency;  Uncoupled  Roll  Channel 
Autopilot;  Classical  Design;  Continuous  Open  Loop 
System;  Circular  Airframe 
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Figure  2.28 


Roll  Angle-Phase  vs  Frequency;  Uncoupled  Roll 
Channel  Autopilot;  Classical  Design;  Continuous 
Open  Loop  System;  Circular  Airframe 
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Figure  2.29 


Roll  Angular  Rate-Gain  vs  Frequency;  Uncoupled  Roll 
Channel  Autopilot;  Classical  Design;  Continuous  Open 
Loop  System  Circular  Airframe 
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Figure  2.30 


Roll  Angular  Rate-Phase  vs  Frequency;  Uncoupled  Roll 
Channel  Autopilot;  Classical  Design;  Continuous  Open 
Loop  System;  Circular  Airframe 
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Figure  2.31 


Roll  Tail  Incidence-Gain  vs  Frequency;  Uncoupled 
Roll  Channel  Autopilot;  Classical  Design; 
Continuous  Open  Loop  System;  Circular  Airframe 
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Figure  2.32 


Roll  Tail  Incidence-Phase  vs  Frequency;  Uncoupled 
Roll  Channel  Autopilot;  Classical  Design; 
Continuous  Open  Loop  System;  Circular  Airframe 
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TABLE  V 

PHASE  CROSSOVER  FREQUENCIES  AND  GAIN  MARGINS;  UNCOUPLED  ROLL 
CHANNEL  AUTOPILOT;  CLASSICAL  DESIGN  CONTINUOUS  OPEN  LOOP  SYSTEM 

CIRCULAR  AIRFRAME 


PHASE  CROSSOVER  GAIN 

FREQUENCY  (rad/sec)       MARGIN  (db) 


ROLL  ANGLE  (ct)  10.0461  60.2165 

ROLL  ANGULAR  RATE  (p)         10.0461  60.2165 

ROLL  TAIL  INCIDENCE  (^  ) 
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4 .   Design  Approach  and  Analysis  of  Discrete  System 

Utilizing  analog-to-digital  conversion  by  the  aid  of 

ORACLS  program  and  for  a  sample  period  of  0.0125  seconds,  a 

seventh-order  discrete  system  of  the  form  x(K  +  1 )=Ax(k)+  Bu(k)  is 

obtained.   The  discrete  plant  system  and  input  matrices  A  and  B 

are  shown  in  Table  VI. 

The  pole-zero  plot  of  Figure  2.33  indicates  that  the 

discrete  open  loop  system  is  also  stable,  since  the  z-plane  poles 

are  : 

-j   =  0.0992805  (roll  angular  rate)  (II. D. 4-1) 

ZX     -    0.835216  +  jO. 309736  (roll  angle)  (II. D. 4-2) 

z5  =  0.83521 6- jO . 309736  (output  of  roll  angle  compensator) 

(II. D. 4-3) 

Z4.    =  0.969087+J0. 0328588  (output  of  rate  compensator  network) 

(II. D. 4-4) 

25  =  0.969087-jO. 0328588  (output  of  pseudo-differentiator) 

(II. D. 4-5) 

£t    -    0.893797  (input  command  in  the  actuator)         (II. D. 4-6) 

Zj    -    0.937075  (roll  tail  incidence)  (II. D. 4-7) 

The  time  response  plots  of  the  roll  angle,  angular  rate 

and  tail  incidence  for  the  discrete  uncoupled  roll  channel  are 

presented  in  Figures  2.34  through  2.36.   A  close  observation  of 

these  plots  indicates  that  they  are  identical  with  those  of  the 

continuous  classical  system  found  in  the  previous  section. 
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Figure  2.33 


Pole-Zero  Plot;  Uncoupled  Roll  Channel;  Classical 
Design;  Discrete  Open  Loop  System;  Circular  Airframe 


81 


o 

0» 

co 

b" 

b" 

<:  ° 

2 

< 

O 

b- 

b" 

CM 

b 

LEGEND 
a    ROLL  ANGLAE.o? 

■4 

© 
©  -f 

0.0 


1.0 


2.0 

TIME  -  SEC 


3.0 


4.0 


Figure    2.3^ 


Roll  Angle  vs  Time;  Uncoupled  Roll  Channel 
Autopilot;  Classical  Design;  Discrete  Open  Loop 
System; 
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Figure    2.35 


Roll  Angular  Rate  vs  Time;  Uncoupled  Roll  Channel 
Autopilot;  Classical  Design;  Discrete  Open  Loop 
System;  Circular  Airframe 
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Figure    2.36 


Roll  Tail  Incidence  vs  Time;  Uncoupled  Roll  Channel 
Autopilot;  Classical  Design;  Discrete  Open  Loop 
System;  Circular  Airframe 


84 


III.   MODERN  CONTROL  DESIGN  AND  ANALYSIS  OF  LINEAR 
UNCOUPLED  LATERAL  AUTOPILOTS 

A.   GENERAL 

The  task  of  this  chapter  is  the  design  and  analysis  of  the 
same  discrete  uncoupled  lateral  channel  autopilots  discussed  in 
the  previous  chapter,  using  different  techniques  which  are  based 
on  modern  control  formulation.   The  difference  in  the  two 
approaches  is  entirely  in  the  design  method  since  the  end  result, 
a  set  of  difference  equations  providing  control,  is  identical. 

Modern  control  theory  is  contrasted  with  the  classical 
control  theory  in  that  the  former  is  applicable  to  multi-input- 
multi-output  systems,  which  may  be  linear  or  nonlinear  time- 
invariant  or  time-varying,  while  the  latter  is  applicable  only  to 
linear  time-invariant  single-input-single-output  systems.   Also, 
modern  control  theory  is  essentially  a  time-domain  approach, 
while  the  conventional  classical  control  theory  is  a  complex 
frequency-domain  approach. 

System  design  in  classical  control  theory  is  based  on  trial- 
and-error  procedures  which,  in  general,  will  not  yield  optimal 
control  systems.   System  design  in  modern  control  theory,  on  the 
other  hand,  enables  the  design  of  optimal  control  systems  of 
great  complexity  and  good  accuracy  with  respect  to  given 
performance  indexes.   In  addition,  design  in  modern  control 
theory  can  be  carried  out  for  a  class  of  inputs  instead  of 
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specific  input  function,  such  as  the  impulse,  step  or  sinusoidal 
functions  and  can  also  include  initial  conditions. 

One  of  the  most  attractive  features  of  modern  control  design 
method  is  that  the  procedure  consists  of  two  independent  steps. 
One  step  assumes  that  all  the  system  states  are  available  for 
feedback  purposes.   In  general,  even  if  this  is  not  a  practical 
enough  assumption  since  it  needs  a  large  number  of  sensors,  it  is 
usually  adopted  in  order  to  accomplish  the  first  design  step, 
namely  the  control-law.   The  remaining  step  is  the  design  of  an 
estimator  which  estimates  the  entire  state  vector,  given 
measurements  of  portion  of  the  state  provided  by  the  system  output 
equation.   The  final  control  algorithm  consists  of  the  control- 
law  and  estimator  combined,  where  the  control-law  calculations 
are  based  on  the  estimated  states  rather  than  the  actual  states. 
This  substitution  is  reasonable  and  the  combined  design  can  give 
closed  loop  characteristics  which  are  unchanged  from  those 
assumed  in  designing  the  control-law  and  estimator  separately. 

B.   DISCRETE  STATE-FEEDBACK  DESIGN 

Considering  the  following  discrete  control  system: 
X(K+1)  =  Ax(k)  +  Bu(k)  (III. B. 1-1) 

Y(k)  =  Hx(k)  (III.B. 1-2) 

the  control-law  design  is  also  referred  to  as  state-feedback 

design  since  it  is  simply  the  feedback  of  a  linear  combination  of 

all  the  system  states,  that  is: 

u(k)  =  Fx(k)  (III.B. 1-3) 

where  F:   control-law  gain  vector 
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Thus,  the  characteristic  equation  of  the  controlled  (closed  loop) 
system  is: 

det  (zI-A+BF)  =  0  (III. B. 1-4) 

The  discrete  state-feedback  design,  providing  that  the  system  is 
controllable,  consists  then  of  finding  the  control-law  gain 
vector  F  so  that  the  roots  of  (III. B.  1-4)  are  in  desirable 
locations . 

A  program  logic  for  computing  the  control-law  gain  vector  via 
the  Ackermann's  formula  is  given  in  Appendix  E  [Ref.  9]. 
Utilizing  this  control  algorithm  a  Fortran  program  was  written 
(Appendix  F)  which  has  as  inputs  the  sample  period,  the  discrete 
plant  system  and  control  input  matrices,  the  s-plane  poles  and 
provides  as  output  the  control-law  gain  vector. 

1 .   Uncoupled  Yaw  Channel  for  Elliptical  Airframe 
a.   Control-Law  Gain  Vector 

Executing  the  Ackermann  Fortran  program  described  in 
Appendix  F  with  inputs: 

(1)  Sample  period  of  0.0125  seconds 

(2)  Discrete  plant  system  and  control  input  matrices  A  and  B 
of  Table  III 

(3)  S-plane  poles  defined  in  equations  (II. C. 3-2)  through 
(II.C.3-8) 

the  following  control-law  gain  vector  for  the  elliptical  airframe 

of  the  uncoupled  yaw  channel  is  obtained: 

F=[-1.0195  -0.109  -0.3492  0.0393  -0.6152  32.4775  -31.0868] 

(III.B. 1-5) 
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b.  Design  Approach  and  Analysis 

The  discrete  state-feedback  designed  yaw  autopilot 
can  be  found  by  introducing  the  control-law  gain  vector  of 
(II. B. 1-5)  into  the  original  system. 

The  pole-zero  plot  of  Figure  3.1  indicates  that  the 

discrete  closed  lop  system  is  stable,  since  the  z-plane  closed 

loop  poles  are: 

Z{      -    0.46461     (yaw   angular    rate)  (III. B. 1-6) 

2^  =  0.898888+ jO. 1 37039  (III. B. 1-7) 

^j  =  0.898888-jO. 137039  (yaw  normal  acceleration)    (III. B. 1-8) 

^4  =  0.998423+jO. 051078  (III. B. 1-9) 

Z<5    -    0.998423-jO. 0510789  (output  of  acceleration  compensator 

network)  (III. B. 1-10) 

Zc    -    0.962974+jO. 0358712  (input  command  in  the  actuator) 

(III.B. 1-11) 

Z>  =  0.962974-jO. 0358712  (yaw  tail  incidence)        (III.B. 1-12) 
The  time  response  plots  of  the  yaw  normal 
acceleration,  angular  rate  and  tail  incidence  are  presented  in 
Figures  3.2  through  3.4.   A  close  observation  of  the  above  pole- 
zero  and  time  response  plots  for  the  discrete  yaw  state-feedback 
design  indicates  that  they  are  identical  with  those  of  the 
discrete  classical  design  found  in  the  previous  chapter. 

c.  Simplified  Design 

The  discrete  state-feedback  designed  autopilot  of  the 
previous  section  can  be  simplified  by  reducing  the  returning  gain 
loops.   This  can  be  accomplished  by  placing  zeros  into  appropriate 
elements  of  the  control-law  gain  vector: 
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F=[-1.0195  -0.109  -0.3492  0  -0.6152  32.4775  -31.0868] 

(III. B. 1-13) 

The  pole-zero  and  time  response  plots  of  the 

resulting  simplified  state-feedback  yaw  autopilot,  shown  in 

Figures  3.5  through  3.8,  do  not  present  significant  differences 

from  the  corresponding  plots  of  the  discrete  classical  design 

apart  from  the  overshoot  which  was  slightly  increased.   The 

discrete  closed  loop  system  is  again  stable,  since  the  z-plane 

closed  loop  poles  are: 

Zi      =  0.467984  )III.B. 1-14) 

Zx     =    0.891328+jO. 126019  (III. B. 1-15) 

z3     =    0.891328-jO. 126019  (III. B. 1-16) 

<*  =  0.998422+jO. 0510788  (III. B. 1-17) 

Z5    =  0.998422-jO. 0510788  (III.B.1.18) 

Zc     -    0.969075  +  jO. 0429695  (III. B. 1-19) 

Zf  =  0.969075-jO. 0429695  (III.  B. 1-20) 

2 .   Uncoupled  Roll  Channel  for  Circular  Airframe 

a.   Control-Law  Gain  Vector 

Following  the  same  procedure  as  in  the  yaw  channel 

case  apart  from  the  use  of  discrete  matrices  A  and  B  from  Table 

III  and  s-plane  poles  from  (II. D. 3-2)  through  (II. D. 3-8),  the 

control-law  gain  vector  for  the  circular  airframe  roll  channel 

was  found  to  be: 

F=[0  0.0004  0  0.0178  -0.0179  -1.499  0.0023]  (III. B. 2-1) 
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Figure  3.1 


Pole-Zero  Plot;  Uncoupled  Yaw  Channel  Autopilot; 
State-Feedback  Design;  Discrete  Closed  Loop  System; 
Elliptical  Airframe 
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Figure    3.2 


Yaw  Normal  Acceleration  vs  Time;  Uncoupled  Yaw 
Channel  Autopilot;  State-Feedback  Design;  Discrete 
Closed  Loop  System;  Elliptical  Airframe 
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Figure    3.3 


Yaw  Angular  Rate  vs  Time;  Uncoupled  Yaw  Channel 
Autopilot;  State-Feedback  Design;  Discrete 
Closed  Loop  System;  Elliptical  Airframe 
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Figure    3.4 


Yaw  Tail  Incidence  vs  Time;  Uncoupled  Yaw  Channel 
Autopilot;  State-Feedback  Design;  Discrete 
Closed  Loop  System;  Elliptical  Airframe 
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Figure  3.5 


Pole-Zero  Plot;  Uncoupled  Yaw  Channel  Autopilot; 
Simplified  State-Feedback  Design;  Discrete  Closed 
Loop  System;  Elliptical  Airframe 
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Figure    3.6 


Yaw  Normal  Acceleration  vs  Time;  Uncoupled  Yaw 
Channel  Autopilot;  Simplified  State-Feedback 
Design;  Discrete  Closed  Loop  System;  Elliptical 
Airframe 
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Figure    3.7 


Yaw  Angular  Rate  vs  Time;  Uncoupled  Yaw 
Channel  Autopilot;  Simplified  State- 
Feedback  Design;  Discrete  Closed  Loop 
System;  Elliptical  Airframe 
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Figure    3-8 


Yaw  Tail  Incidence  vs  Time;  Uncoupled  Yaw 
Channel  Autopilot;  Simplified  State- 
Feedback  Design;  Discrete  Closed  Loop 
System;  Elliptical  Airframe 
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b.  Design  Approach  and  Analysis 

The  discrete  state-feedback  designed  roll  autopilot 
can  be  found  by  introducing  the  control-law  gain  vector  of 
(III. B. 2-1)  into  the  original  system. 

The  pole-zero  plot  of  Figure  3-9  indicates  that  the 

discrete  closed  loop  system  is  stable,  since  the  z-plane  closed 

loop  poles  are: 

2,   =  0.0989714  (roll  angular  rate)  (III. B. 2-2) 

ZZ      =  0.832909+jO. 315088  (roll  angle)  (III. B. 2-3) 

^5  =  0.832909-jO. 315088  (output  of  roll  angle  compensator) 

(III. B. 2-4) 

Z*.     =  0 . 969236+ jO. 0333773  (output  of  rate  compensator  network) 

(III. B. 2-5) 

Zs    =  0.969236-jO. 0333773  (output  of  pseudo-differentiator) 

(III. B. 2-6) 

Z£     =  0.896583  (input  command  in  the  actuator)        (III. B. 2-7) 

Zf    =  0.938233  (roll  tail  incidence)  (III. B. 2-8) 

The  time  response  plots  of  the  roll  angle,  angular 

rate  and  tail  incidence  are  presented  in  Figures  3.10  through 

3.12.   A  close  observation  of  the  above  pole-zero  and  time 

response  plots  for  the  discrete  roll  state-feedback  design 

indicates  that  they  are  identical  with  those  of  the  discrete 

classical  design  found  in  the  previous  chapter. 

c.  Simplified  Design 

The  state-feedback  roll  designed  autopilot  can  be 
simplified  by  reducing  the  returning  gain  loops  as  follows: 
F=[0  0  0  0.178  -0.0179  -1.499  0]  (III. B. 2-9) 
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The  pole-zero  and  time  response  plots  of  the 
resulting  simplified  state-feedback  roll  autopilot,  shown  in 
Figures  3.13  through  3.16,  do  not  present  significant  differences 
from  the  corresponding  plots  of  the  discrete  classical  design. 
The  discrete  closed  loop  system  is  again  stable,  since  the  z- 
plane  closed  loop  poles  are: 

/,   =  0.0992745  (III. B. 2-10) 

*1  =  0.836632+jO. 305719  (III. B. 2-11) 

Z5  =  0.836632-jO. 305719  (III. B. 2-12) 

^4  =  0.967647+jO. 0334948  (III. B. 2-13) 

ZS   =  0.967647-jO. 0334948  (III. B. 2-14) 

Z€    =  0.895937  (III. B. 2-15) 

^7  =  0.938256  (III. B. 2-16) 

C.   DISCRETE  ESTIMATOR  DESIGN 

The  state-feedback  design  discussed  in  the  last  section 
assumed  that  all  system  states  were  available  for  feedback 
purposes.   Since  the  state  vector  is  not  always  accessible  to 
direct  measurement,  an  estimator  is  going  to  be  introduced  in 
this  section  as  an  additional  dynamic  design  in  order  to 
implement  control  to  the  original  system.   The  estimator  design 
method  consists  mainly  of  determining  algorithms  which  will 
reconstruct  all  the  states,  given  measurements  of  a  portion  of 
them. 
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Figure  3.9 


Pole-Zero  Plot;  Uncoupled  Roll  Channel  Autopilot; 
State-Feedback  Design;  Discrete  Closed  Loop  System; 
Circular  Airframe 
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Figure    3.10 


Roll  Angle  vs  Time;  Uncoupled  Roll  Channel  Autopilot; 
State-Feedback  Design;  Discrete  Closed  Loop  System; 
Circular  Airframe 
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Figure    3. 1 1 


Roll  Angular  Rate  vs  Time;  Uncoupled  Roll  Channel 
Autopilot;  State-Feedback  Design;  Discrete  Closed 
Loop  System;  Circular  Airframe 
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Figure    3.12 


Roll  Tail  Incidence  vs  Time;  Uncoupled  Roll  Channel 
Autopilot;  State-Feedback  Design;  Discrete  Closed 
Loop  System;  Circular  Airframe 
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Figure  3.13 


Pole-Zero  Plot;  Uncoupled  Roll  Channel  Autopilot; 
Simplified  State-Feedback  Design;  Discrete  Closed 
Loop  System;  Circular  Airframe 
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Figure    3.14 


Roll  Angle  vs  Time;  Uncoupled  Roll  Channel  Autopilot; 
Simplified  State-Feedback  Design;  Discrete  Closed 
Loop  System;  Circular  Airframe 
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Figure  3.15 


Roll  Angular  Rate  vs  Time;  Uncoupled  Roll  Channel 
Autopilot;  Simplified  State-Feedback  Design; 
Discrete  Closed  Loop  System;  Circular  Airframe 


106 


0.0 


1.0 


2.0 

TIME  -  SEC 


3.0 


4.0 


Figure    3. 1 6 


Roll  Tail  Incidence  vs  Time;  Uncoupled  Roll  Channel 
Autopilot;  Simplified  State-Feedback  Design; 
Discrete  Closed  Loop  System;  Circular  Airframe 
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Considering  the  same  as  in  the  state-feedback  case  discrete 
control  system,  a  prediction  estimator  defined  by  the  following 


equation  is  introduced: 

x(K+1)  =  Ax(k)  +  Bu(k)  +  K[X(k)  -  Hx(k)] 

A 

where  x:  estimate  state  vector 


(III.C. 1-1  ) 


K:  estimate  gain  vector 
In  this  closed  loop  estimator,  shown  in  Figure  3.17,  the 
difference  between  the  measured  and  estimated  output  is  fed  back 
and  the  model  is  constantly  corrected  with  this  error  signal 
which  is  defined  as  x=x-x.   The  difference  equation  describing 
the  behavior  of  the  error  is  obtained  by  subtracting  equation 
from  the  actual  plant  output  equation  (II. B. 1-2): 

x(K+1)  =  [A  -  KH]  x(k)  (III.C. 1-2) 

Thus  the  characteristic  equation  of  the  controlled  (closed  loop) 
system  is: 

det  (zl  -  A  +  KH)  =  0  (III.C. 1-3) 

The  discrete  estimator  design,  providing  that  the  system  is 
observable,  consists  then  of  finding  the  estimator  gain  vector  K 
so  that  the  roots  of  (III.C. 1-3)  are  at  desirable  locations. 

The  estimator  gain  vector  can  be  obtained  again,  as  in  the 
case  of  the  state-feedback  design,  by  application  of  the 
Ackermann  program  of  Appendix  F  with  inputs  the  sample  period, 
the  transposes  of  the  discrete  plant  and  output  matrices,  and 
faster  s-plane  poles  from  those  of  the  continuous  system. 
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1 .   Uncoupled  Yaw  Channel  for  Elliptical  Airframe 
a.   Estimator  Gain  Vector 

Executing  the  Ackermann  Fortran  program  (Appendix  F) 
with  inputs: 

(1)  Sample  period  of  0.0125  seconds 

(2)  The  transposes  of  the  discrete  yaw  system  plant  and  output 

T         T 

matrices,  that  is  A1  and  H-i  ,  where: 


i-J. 


1 

0 
0 
0 

0 
0 
0 


(III.C. 1-4) 


(3)   S-plane  poles  slightly  faster  than  those  of  the  continuous 
open  loop  system,  that  is  with  more  negative  real  parts: 


S,  =  -174.386 

5^  =  -6. 1 142+J10.6396  • 

S<  =  -6. 1 1 42- j 1 0 . 6396 

i+  =  -0.0227475+J4. 08919 

S5  =  -0.0227475-J4. 08919 

Sg  =  -2.9396+J2.99929 

$r  =  -2.9396-J2.99929 


(III.C. 1-5) 
(III.C. 1-6) 
(III.C. 1-7) 
(III.C. 1-8) 
(III.C. 1-9) 
(III.C. 1-10) 
(III.C. 1-1 1 ) 


the  transpose  of  the  estimator  gain  vector  for  the  elliptical 
airframe  of  the  uncoupled  yaw  channel  is  calculated  as  output, 
from  which  the  following  K  is  obtained: 
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K  = 


-0.003^ 
0.0008 

-0.0003 
-0.0177 
0.0005 
-0.0017 
-0.0240 


(III.C. 1-12) 


b.   Design  Approach  and  Analysis 

The  discrete  estimator  designed  yaw  autopilot  can  be 
found  by  introducing  the  estimator  gain  vector  of  (III.C. 1-12) 
into  the  original  system. 

The  pole-zero  plot  of  Figure  3.18  indicates  that  the 
discrete  closed  loop  system  is  stable,  since  the  z-plane  closed 
loop  poles  were  found  to  be: 
Z,   =  0.117115  (yaw  angular  rate) 
?2     -    0.918174+jO. 122926 

^5  =  0.9181 74- JO . 122926  (yaw  normal  acceleration) 
zA    =  0.998399+jO. 0510958 


(III.C. 1-13) 
(III.C. 1-14) 
(III.C. 1-15) 
(III.C. 1-16) 


£<:    =    0.998399-jO. 0510958    (output    of    acceleration    compensator 


network) 


(III.C. 1-17) 


Zc      -    0.963249+jO .0361 306  (input  command  in  the  actuator) 

(III.C. 1-18) 

Zj    =  0.963249-jO. 0361306  (yaw  tail  incidence)        (III.C. 1-19) 

The  time  response  plots  of  the  yaw  normal 
acceleration,  angular  rate  and  tail  incidence  are  presented  in 
Figures  3.19  through  3.21.   A  close  observation  of  the  above  pole- 
zero  and  time  response  plots  of  the  discrete  yaw  estimator  design 
indicates  that  they  are  very  close  to  those  of  the  discrete 
classical  design  found  in  the  previous  chapter. 
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Figure    3.18 


Pole-Zero  Plot;  Uncoupled  Yaw  Channel  Autopilot; 
Estimator  Design;  Discrete  Closed  Loop  System; 
Elliptical  Airframe 
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c.   Simplified  Design 

The  discrete  estimator  designed  autopilot  of  the 
previous  section  can  be  simplified  by  reducing  the  returning  gain 
loops.   This  can  be  accomplished  by  placing  zeros  into 
appropriate  elements  of  the  estimator  gain  vector: 


K  = 


-0.0034 

0 

0 
-0.0177 

0 
-0.0017 
-0.0240 


(III.C. 1-20) 


The  pole-zero  and  time  response  plots  of  the 
resulting  simplified  estimator  yaw  autopilot,  shown  in  Figures 
3.22  through  3.25,  do  not  present  significant  differences  from 
the  corresponding  plots  of  the  discrete  classical  design,  apart 
from  the  overshoot  which  was  slightly  increased.   The  discrete 
closed  loop  system  is  again  stable,  since  the  z-plane  closed  loop 
poles  are: 

Z,   =0.117117  (III.C. 1-21  ) 

ZA     =  0.918434+jO. 122898  •  (III.C. 1-22) 

Z2    -    0.918434-jO. 122898  (III.C. 1-23) 

z^  =  0.998422+J0. 0510954  (III.C. 1-24) 

^5  =  0.998422-jO. 0510954  (III.C. 1-25) 

ze    =  0.962965+jO. 0357874  (III.C. 1-26) 

Z7   =  0.962965-jO. 0357874  (III.C. 1-27) 
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2 .   Uncoupled  Roll  Channel  for  Circular  Airframe 
a.   Estimator  Gain  Vector 

Following  the  same  procedure  as  in  the  yaw  channel 
case  apart  from  the  use  of: 

(1)  The  transpose  of  the  discrete  roll  system  plant  matrix 

(2)  S-Plane  poles  slightly  faster  than  those  of  the  continuous 
open  loop  roll  system,  that  is: 

5,   =  -184.795  (III. C. 2-1 ) 

5Z     -    -9.26097  +  J28.4098  (III. C. 2-2) 

%  =  -9.26097-j28.4098  (III. C. 2-3) 

%  =  -2.47608+J2.71152  (III. C. 2-4) 

S5  =  -2 . 47608- j2 .71 1 52  (III. C. 2-5) 

$c    =  -8.99209  (III. C. 2-6) 

S7    =  -5.20032  (III. C. 2-7) 

the  estimator  gain  vector  for  the  circular  airframe  roll  channel 

is  obtained: 


0.004 
0 
0 

-0.0002 
0.0003 
0 
-0.0492 


(III. C. 2-8) 


b.   Design  Approach  and  Analysis 

The  discrete  estimator  designed  roll  autopilot  can  be 
found  by  introducing  the  estimator  gain  vector  of  (III. C. 2-8) 
into  the  original  system. 
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Figure    3.19 


Yaw  Normal  Acceleration  vs  Time;  Uncoupled  Yaw 
Channel  Autopilot;  Estimator  Design;  Discrete 
Closed  Loop  System;  Elliptical  Airframe 
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Figure    3.20 


Yaw  Angular  Rate  vs  Time;  Uncoupled  Yaw 
Channel  Autopilot;  Estimator  Design; 
Discrete  Closed  Loop  System;  Elliptical 
Airframe 
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Figure    3.21 


Yaw  Tail  Incidence  vs  Time;  Uncoupled  Yaw 
Channel  Autopilot;  Estimator  Design; 
Discrete  Closed  Loop  System;  Elliptical 
Airframe 
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Figure  3.22 


Pole-Zero  Plot;  Uncoupled  Yaw  Channel 
Autopilot;  Simplified  Estimator  Design; 
Discrete  Closed  Loop  System;  Elliptical 
Airframe 
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Figure    3.23 


Yaw  Normal  Acceleration  vs  Time;  Uncoupled  Yaw 
Channel  Autopilot;  Simplified  Estimator  Design; 
Discrete  Closed  Loop  System;  Elliptical  Airframe 
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Figure    3.24 


Yaw  Angular  Rate  vs  Time;  Uncoupled  Yaw 
Channel  Autopilot;  Simplified  Estimator 
Design;  Discrete  Closed  Loop  System; 
Elliptical  Airframe 
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Figure    3.25 


Yaw  Tail  Incidence  vs  Time;  Uncoupled  Yaw 
Channel  Autopilot;  Simplified  Estimator 
Design;  Discrete  Closed  Loop  System; 
Elliptical  Airframe 
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The  pole-zero  plot  of  Figure  3.26  indicates  that  the 
discrete  closed  loop  system  is  stable,  since  the  z-plane  closed 
loop  poles  were  found  to  be: 


Z,  =  0.0967149 

^  =  0.830876+jO. 31913 

z5  =  0 .  830876- jO .31913 

Z4  =  0.970259+jO. 0325175 

Zs    =  0.970259-jO. 0325179 

Zc    =  0.894774 

l~l   =    0.937225 


(III. C. 2-9) 
(III. C. 2-10) 
(III. C. 2-1 1  ) 

(III. C. 2-13) 
(III. C. 2-14) 
(III. C. 2-15) 
(III. C. 2-16) 


The  time  response  plots  of  the  roll  angle,  angular 
rate  and  tail  incidence  are  presented  in  Figures  3.27  through 
3.29.   A  close  observation  of  the  above  pole-zero  and  time 
response  plots  for  the  discrete  roll  estimator  design  indicates 
that  they  are  very  close  to  those  of  the  discrete  classical  design 
found  in  the  previous  chapter. 
c.   Simplified  Design 

The  estimator  roll  designed  autopilot  can  be 
simplified  by  reducing  the  returning  gain  loops  as  follows: 


K  = 


0 
0 
0 

0.0178 

-0.0179 
-1 .499 
0 


(III. C. 2-17) 


The  pole  and  time  response  plots  of  the  resulting 
simplified  estimator  roll  autopilot,  shown  in  Figures  3-30 
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through  3-33,  do  not  present  significant  differences  from  the 
corresponding  plots  of  the  discrete  classical  design.   The 
discrete  closed  loop  system  is  again  stable,  since  the  z-plane 
closed  loop  poles  are: 

H,  =  0.0992745  (III. C. 2-18) 

^Z    =  0.836632+jO. 305719  (III. C. 2-19) 

Zj     =  0 . 836632- jO . 30571 9  (III. C. 2-20) 

£4  =  0.967647+jO. 033^948  (III. C. 2-21) 

%  =  0.967647-jO. 0334948  (III. C. 2-22) 

*c     =  0.895937  (III. C. 2-23) 

^7  =  0.938256  (III. C. 2-24) 
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Figure  3.26 


Pole-Zero  Plot;  Uncoupled  Roll  Channel  Autopilot; 
Estimator  Design;  Discrete  Closed  Loop  System; 
Circular  Airframe 
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Figure    3-27 


Roll  Angle  vs  Time;  Uncoupled  Roll  Channel  Autopilot; 
Estimator  Design;  Discrete  Closed  Loop  System; 
Circular  Airframe 
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Figure    3.28 


Roll  Angular  vs  Time;  Uncoupled  Roll  Channel 
Autopilot;  Estimator  Design;  Discrete  Closed  Loop 
System;  Circular  Airframe 


126 


LEGEND 
D    ROLL  TAIL  INCIDENCE. (5R 


o.o 


1.0 


2.0 

TIME  -  SEC 


3.0 


4.0 


Figure    3.29 


Roll  Tail  Incidence  vs  Time;  Uncoupled  Roll  Channel 
Autopilot;  Estimator  Design;  Discrete  Closed  Loop 
System;  Circular  Airframe 


127 


o 

o 

(9 

o  - 

b 

N 

X 

O 

v.    / 

.  .x 

c. 

o 

•      X 

;vx 

1 

. 

X 

1 

• 

1 

09 

1 

o 
1 

1 

1 

-1.00 


-0.75 


-0  50 


-0.25 


0.00 


0.25 


0.50 


0.75 


1.00 


Figure  3-30 


Pole-Zero  Plot;  Uncoupled  Roll  Channel  Autopilot; 
Simplified  Estimator  Design;  Discrete  Closed  Loop 
System;  Circular  Airframe 
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Figure    3-31 


Roll  Angle  vs  Time;  Uncoupled  Roll  Channel  Autopilot; 
Simplified  Estimator  Design;  Discrete  Closed  Loop 
System;  Circular  Airframe 
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Figure  3.32 


Roll  Angular  Rate  vs  Time;  Uncoupled  Roll  Channel 
Autopilot;  Simplified  Estimator  Design;  Discrete 
Closed  Loop  System;  Circular  Airframe 
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Figure    3.33 


Roll  Tail  Incidence  vs  Time;  Uncoupled  Roll  Channel 
Autopilot;  Simplified  Estimator  Design;  Discrete 
Closed  Loop  System;  Circular  Airframe 
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IV.   MODERN  CONTROL  DESIGN   AND  ROBUSTNESS  ANALYSIS  OF 
COUPLED  PITCH  AND  ROLL  CHANNEL  AUTOPILOT 
CIRCULAR  AIRFRAME 

A.   GENERAL 

The  present  chapter  deals  with  the  modern  control  design  and 
robustness  analysis  of  the  discrete  coupled  pitch  and  roll 
channel  autopilot  for  the  circular  airframe  configuration.   The 
continuous  open  loop  coupled  autopilot  whose  plant  system  and 
input  matrices  are  presented  in  Appendix  G  is  obtained  by 
coupling  the  linear  uncoupled  pitch  [Ref.  10]  and  roll  (Table  IV) 
channels.   Then,  utilizing  analog-to-digital  conversion  by  the 
aid  of  ORACLS  program  and  for  a  sample  period  of  0.0125  seconds, 
the  seventeenth-order  discrete  coupled  system  with  matrices  shown 
in  Appendix  H  is  formulated.   Next,  introducing  the  control-law 
and  estimator  designs  were  obtained  and  analyzed  in  terms  of  their 
transient  responses  and  the  application  of  the  POPLAR  design 
program  [Ref.  7].   The  POPLAR  program  is  applied  in  order  to 
employ  singular  value  analysis  and  the  use  of  an  optimization 
routine  to  aid  in  pole  placement  control  design  of  the  above 
discussed  linear  multivariable  systems.   The  robustness  of  the 
system  is  also  considered  by  establishing  singular  value  levels 
which  correspond  to  multiloop  gain  and  phase  margins  determined 
from  the  universal  gain  phase  system  diagram  developed  by  Newsom 
and  Mukhapadhyay  at  NASA  Langley. 
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B.   DISCRETE  COUPLED  STATE-FEEDBACK  DESIGN 


1 .   Design  Approach  and  Analysis 

The  discrete  coupled  state-feedback  designed  autopilot  is 
formulated  by  introducing  into  the  original  control  system  of 
Appendix  H  the  following  combined  pitch  [Ref.  10]  and  roll 
(III. B. 1-5)  control-law  gain  vector. 


O.OI88 

-o.cw 

-  O.CC4-Z 
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O.OOOi 

\.rrsz 
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O.  OCo4 

o 

O.OIT9 

-CCI79 

-14-99 

C. 

0013 

(IV. B. 1-1 ) 
The  pole-zero  plot  of  Figure  4.1  indicates  that  the 
discrete  coupled  system  is  marginally  stable,  since  the  z-plane 
closed  loop  poles  are: 

£,  =  0.  13446  +  jO. 0320974  (IV. B. 1-2) 

^  =  0. 1 3446- jO. 0320974  (IV. B. 1-3) 

^3  =  0.89699+jO. 0906457  (IV. B. 1-4) 

Z+     =  0.89699- j0. 0906457  (IV. B. 1-5) 

%  =  0.955624+jO. 0299994  (IV. B. 1-6) 

£c     =  0.953624-jO. 0299994  (IV. B. 1-7) 

^7  =  0.998202  (IV. B. 1-8) 

(IV. B. 1-9) 
(IV. B. 1-10) 


**  =  1 

*V  =  0.995608+jO. 0830713 
^10   =  0.995608-J0.0830713 


(IV. B. 1-11) 
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Figure  4  .  1 


Pole-Zero  Plot;  Coupled  Pitch  and  Roll  Channel 
Autopilot;  State-Feedback  Design;  Discrete 
Closed  Loop  System;  Circular  Airframe 
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2,,    =    0.0993086  (IV. B.  1-12) 

Ba   =    0.836555+jO. 306107  (IV. B. 1-13) 

2|3  =    0.836553-jO. 306107  (IV. B. 1-14) 

£,4.s  0.967788  +  jO. 033^452  (IV. B. 1-15) 

2,5  =  0.967788-jO. 0334452  (IV. B. 1-16) 

2,fi  =  0.89573+jO. 938133  (IV. B. 1-17) 

2,7  =  0.89573-jO. 938133  (IV. B. 1-18) 

The  time  response  plots  of  the  roll  angle,  angular  rate 
and  tail  incidence  are  presented  in  Figures  4.2  through  4.4.   A 
close  observation  of  these  plots  indicate  that  they  are  identical 
with  those  of  the  discrete  uncoupled  state-feedback  design  of  the 
previous  chapter. 

2 .   Robustness  Analysis 

Executing  the  POPLAR  design  program  of  [Ref.  7]  with 
inputs  the  data  presented  in  Appendix  I,  the  minimum  additive 
input  (MIN  ADD  IN  SV)  and  output  (SVADMO)  singular  values  were 
computed  from  a  frequency  range  from  0  to  200  rad/sec. 

Figures  4.5  and  4.6  which  are  plots  of  SVADMO  and  MIN  ADD 
IN  SV  versus  frequency  indicate  that  the  discrete  coupled  state- 
feedback  design  is  robust.   It  is  noted  that  for  very  low 
frequencies  the  values  of  SVADMO  are  above  0.8680. 

Finally,  in  terms  of  optimization  results,  the  ordered 
computed  eigenvalues  are: 

2,   =  0.07717  (IV. B. 2-1 ) 

Zz     -    0.09855  (IV. B. 2-2) 

Z5  =  0.25923  (IV. B. 2-3) 
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Figure    4.2 


Roll  Angle  vs  Time;  Coupled  Pitch  and  Roll  Channel 
Autopilot;  State-Feedback  Design;  Discrete  Closed 
Loop  System;  Circular  Airframe 
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Figure    4.3 


Roll  Angular  Rate  vs  Time;  Coupled  Pitch  and  Roll 
Channel  Autopilot;  State-Feedback  Design;  Discrete 
Closed  Loop  System;  Circular  Airframe 
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Figure  4.4 


Roll  Tail  Incidence  vs  Time;  Coupled  Pitch  and  Roll 
Channel  Autopilot;  State-Feedback  Design;  Discrete 
Closed  Loop  System;  Circular  Airframe 
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Figure    4.5 


SVADMO  vs  Frequency;  Coupled  Pitch  and  Roll 
Channel  Autopilot;  State-Feedback  Design; 
Discrete  Closed  Loop  System;  Circular  Airframe 
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MIN  ADD  IN  SV  vs  Frequency;  Coupled  Pitch  and  Roll 
Channel  Autopilot;  State-Feedback  Design;  Discrete 
Closed  Loop  System;  Circular  Airframe 
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ftA     =  0.069508  (IV. B. 2-4) 

£5  =  0  .  83707- j0  .3014  (IV. B. 2-5) 

■2£  =  0.83707  +  J0.3014  (IV. B. 2-6) 

27  =  0.89844  (IV. B. 2-7) 

Ls    =  0.90064  (IV. B. 2-8) 

±c,    =  0.43996  (IV. B. 2-9) 

2,0  =  0.46207  (IV. B. 2-10) 

2n  =  0.96562-jO. 03401  (IV. B. 2-11) 

Za   =  0.96502+jO. 03401  (IV. B. 2-12) 

2l5  =  0.99559-jO. 08307  (IV. B. 2-13) 

2I4.  =  0.99559  +  jO. 08307  (IV. B. 2-14) 

2l5  =  0.99820  (IV. B. 2-15) 

£IC   =  1  (IV. B. 2-16) 

2f7=  1.06546  (IV. B. 2-17) 

C.   DISCRETE  COUPLED  ESTIMATOR  DESIGN 

1 .   Design  Approach  and  Analysis 

Following  the  same  procedure  as  in  the  case  of  the  state- 
feedback  but  for  the  coupled  estimator  gain  vector  the  pole-zero 
and  time  response  plots  of  the  coupled  estimator  design  are 
obtained.   The  time  responses  are  again  identical  with  those  of 
the  discrete  uncoupled  estimator  design. 

2  .   Robustness  Analysis 

Executing  the  POPLAR  design  program  of  [Ref.  7],  Figures 
4.7  and  4.8  are  obtained  which  prove  the  robustness  of  the 
system.   It  is  noted  that  for  very  low  frequencies  the  values  of 
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SVADMO  are  above  0.83976.   Finally,  in  terms  of  optimization 
results  the  ordered  computed  eigenvalues  are  the  same  as  in  the 
coupled  state-feedback  design. 
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Figure  4.7 


SVADMO  vs  Frequency;  Coupled  Pitch  and  Roll 
Channel  Autopilot;  Estimator  Design;  Discrete 
Closed  Loop  System;  Circular  Airframe 
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MIN  ADD  IN  SV  vs  Frequency;  Coupled  Pitch  and  Roll 
Channel  Autopilot;  Estimator  Design;  Discrete 
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V.   CONCLUSIONS  AND  RECOMMENDATIONS 

A.  CONCLUSIONS 

The  goal  of  the  present  thesis  was  the  design  and  analysis  of 
discrete  lateral  autopilots  for  application  to  BTT  missiles.  The 
following  are  the  principal  conclusions  based  on  this  work. 

1.  The  continuous  and  discrete  classical  designed  autopilots 
were  proved  to  have  identical  performances  for  the  two 
lateral  channels. 

2.  The  state-feedback  and  estimator  autopilots  were  introduced 
as  additional  dynamic  designs  in  order  to  implement  control 
to  the  original  system.   Both  designs,  analyzed  in  terms  of 
their  transient  responses,  were  found  to  meet  the  desired 
requirements . 

3.  The  simplified  state-feedback  and  estimator  designs  reduced 
some  of  the  returning  gain  loops,  making  the  system 
simpler,  without  any  significant  effects  on  the  system's 
performance . 

4.  The  performance  of  the  coupled  pitch  and  roll  channel 
autopilot  was  found  to  be  satisfactory  and  the  overall 
system  proved  to  be  robust. 

B.  RECOMMENDATIONS 

In  order  to  improve  the  simplicity  of  the  overall  system, 
more  returning  gain  loops  of  the  state-feedback  and  estimator 
design  must  be  eliminated.   A  further  investigation  then  must  be 
conducted  in  order  to  examine  if  the  performance  of  the  resulting 
design  remains  unchanged. 
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APPENDIX  A 
DESIGN  REQUIREMENTS  FOR  UNCOUPLED  AUTOPILOTS 

The  requirements  for  the  classical  design  method  of  the 
uncoupled  channel  autopilots  [Ref.  6]  are  the  following: 

1 .  High  Frequency  Attentuation  in  Actuator  Command  Branch 

a.  Uncoupled  Yaw  Channel 

It  must  be  >y   15  db  at  100  rad/sec  and  zero  angle-of- 
attack  and  sideslip.   This  requirement  will  provide  sufficient 
high  frequency  attenuation  for  >y   30  Hz  actuator  and  for  body 
bending  modes  when  high  frequency  filters  are  added,  but  it 
limits  the  ability  of  the  yaw  autopilot  to  minimize  sideslip 
angle  . 

b.  Uncoupled  Roll  Channel 

It  must  be  >y   15  db  at  100  rad/sec  and  zero  angle-of- 
attack.   This  will  provide  sufficient  high  frequency  attenuation 
for  >,  30  Hz  actuator  and  for  elastic  modes  when  high  frequency 
filters  are  added,  but  this  requirement  limits  the  speed  of  roll 
angle  response. 

2 .  Relative  Stability  for  Both  Lateral  Channels 

Gain  margins  >6  db,  phase  margins  ^30°  with  a  goal  of 
12  db  and  50°. 

3  .   Acceleration  Time  Response 
a.   Uncoupled  Yaw  Channel 

(1)  63$  time  constant  of  approximately  0.4  seconds. 
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(2)  Overshoot  <<:  10%. 

(3)  Steady-state  error  need  not  be  zero 
Uncoupled  Roll  Channel 

(1)  63%  time  constant  of  0.5  seconds. 

(2)  Overshoot  <:  10%. 

(3)  Zero  steady-state  roll  angle  error. 
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APPENDIX  B 
AERODYNAMIC  DATA 

The  overall  classical  design  developed  by  Arrow  [Ref.  6]  was 
performed  for  the  selected  flight  condition  of  Mach  number  3-95 
at  60000  feet  altitude.   The  corresponding  aerodynamic  data 
presented  below  were  taken  or  derived  from  the  ICAO  standard 
atmosphere  tables. 

TABLE  VII 
AERODYNAMIC  DATA  (M=3.95,  H=60kft) 

Temperature,  T  (°R)  389.988 

Sonic  Velocity,  a  (ft/sec)  968.47 

Pressure,  p  (lb/ft2)  149.78 

Density,  (lb-sec2/ft)  0.0002238 

Velocity,  V  (ft/sec)  3825.4565 

Dynamic  Pressure,  q  (lb/ft2)  1637.145 
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APPENDIX  C 
LINEARIZED  AERODYNAMIC  DERIVATIVES 

The  linearized  aerodynamic  derivatives  at  the  selected  flight 
condition  and  about  a  zero  trim  angle-of-attack  are  provided 
below  for  both  airframe  configurations  [Ref.  6]. 

TABLE  VIII 

LINEARIZED  AERODYNAMIC  DERIVATIVES 
(M=3.95,  H=60kft;  a=0°) 

CIRCULAR  ELLIPTICAL 

cv  -0.065  -0.043 


nP 


'o 


-0.025  0.024 

0.021  0.016 

Y  -0.050  -0.042 


u£r  0.031  0.023 

"ft 
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APPENDIX  D 
MISSILE  SIZING  AND  MASS  PROPERTIES 

In  order  to  provide  a  realistic  missile  based  on  the 
aerodynmamically  tested  configuration  concepts  [Ref.  5],  the 
models  were  assumed  to  be  1/6-scale  and  the  mass  properties  were 
developed  corresponding  to  mass  distribution  which  might  be 
expected  for  missiles  of  this  size.   All  the  geometric  and  mass 
properties  are  presented  in  the  following  table. 

TABLE  IX 
GEOMETRIC  AND  MASS  PROPERTIES  OF  MISSILE  CONFIGURATIONS 


Length,  1  (in) 

Max.  Diameter  (in) 

Max.  Major  Axis  (in) 

Max.  Minor  Axis  (in) 

e.g.  distance  from  L.E.  (in) 

Reference  Length,  d  (ft) 

Reference  Area,  S  (ft) 

Weight,  W  (lb) 

Ixx  (slug-ft2) 

Izz  (slug-ft2) 


Circular 


168 


24 


Elliptical 
168 


41  .57 

13-86 

100.8(0, 

.6  ) 

100.8(0.6  ) 

2 

2 

TT 

TT 

2525 

2475 

40 

40 

810 

853 
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APPENDIX  E 
PROGRAM  LOGIC  FOR  APPLICATION  OF  ACKERMANN'S  FORMULA 

The  program  logic  for  computing  the  control-law  gain  vector  F 

T 

or  the  transpose  of  the  estimator  gain  vector  K1 ,  taken  from 
[Ref.  9],  is  given  in  the  following  table. 

TABLE 
PROGRAM  LOGIC  FOR  APPLICATION  OF  ACKERMANN'S  FORMULA 


I  Read  in  <\k  V.  T.  and  ,Y„.  the  number  of  slates. 

2.  Comment:  first  we  will  read  in  the  desired  pole  locations  in  the  \-plane.  convert  them 
to  ;-plane  polynomial  coefficients,  and  construct  ui<t>). 

3.  I  —  identity  matrix.  \t  x    \ 
4  ALPHA  —  I 

5.  k  —  1 

b.  If  k     ■   V, .  go  to  step  18. 

7.  Read  in  pole  location  k  as  u  -  jh. 

8.  If  h  =  0.  yo  to  step  14. 

9.  A,  > 2  e\p<</7  >  cos  hT 

10.  A,  —  e\pt2,/7  i 

II  ALPHA  —ALPHA   -   i'I>       <1>  -    \  ,<!»  -  ,\.\) 

12.  k  —  k   ~  2 

13.  Go  to  step  6 

14.  A,  «—  e\pi(;7  ) 

15.  ALPHA  —  ALPHA   >   («l>  -  A,   >    I) 

16.  A.  *-  k  ->   1 

17.  Go  to  step  6. 

IK.  Comment:  nou  we  construct  the  controllability  matrix. 

19.  C  —  I 

20.  E  —  r 

21.  A  —  I 

22.  If  k  >  A',,  go  to  step  28. 

23.  Comment:  replace  column  k  of  C  bv  E. 

24.  C[      :  k]  —  E 

25.  A  —  A  -   I 

26.  E  —  <l>   •    E 

27.  Go  to  step  22. 

28.  Comment:  mm  solve  for  the  control  law.  hist  form  c7„  as  the  last  row  of  I. 

29.  E  —  i[.V„:     ] 

30.  Solve  KC  =  E  for  B. 

31.  K  =  B   ■    ALPHA 

32.  END 
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APPENDIX  F 


ACKERMANN  FORTRAN  PROGRAM 


FILE:     ACKERMAN    /.ATrlV 


Al 


sjca 


c 
c 

c 


INTEGER  NS, * , I  A ,  IER,  I  ,J,NI 

REAL  Ph( 20,20 ) , GA I  2  0,2  J  J. A  It  20, 20), ALPHA (20,20)  ,*M  40) 

REAL  A2        ,t  (23,20)  ,01  20,2  J)  ,3(20)  ,  BF  U  J, 20)  ,  JK(20,2  J) 

REAL       A,3,T,     R,  ATEMk  I  20,  2  0) 

REAL  «1  , TG12U, 20) , Til  20.20) ,TSl 20,20) ,u(20, 20), TM(20, 20) 


INPUT    OF    NS=    NUMBER    OF    STATES 
T=    SaMPLc    TIME 

1     ENTER    NUMBER    OF    STATES.' 
IS 
•    ENTER    SAMPLE    TIME    ' 


IN 


103 

101 
100 

IN 


105 

200 
BU 
At, 


3  00 


C 

c 
c 


16 


PRIN 
REAl) 
PRIN 
READ 

PUT  J 
00  1 
DC  I 
WRI  T 
FORM 
READ 
CONT 
CGNT 

PUT  0 
DG  2 
J  = 
WRI  T 
FORM 
REAO 
CUNT 

1LDIN 

PHA=C 
DO  1 
DG  5 
IF  ( 
All  I 
ALPri 
C(I  , 
GC  T 
All  I 
ALPH 
C(l  , 
CCNT 
CCNT 
K.=  l 
IF     ( 


T 

,    N 
T, 
,    T 
F     P 
00 
01 
E(  6 
ATI 

,?H 
1NU 
INU 
F  G 
00 
1 

E(6 
AT( 
,GA 
INU 
G  T 
=  i 
1  = 
J  = 
I  .E 
,J) 
A(  I 
J)  = 
Li  5 
,J) 
Al  1 
J)  = 
INU 
INU 


HI  =  A     MATRIX 
1=1, NS 
J=L,NS 
,103)     I, J 

<,X,  '  ENTEk    PHI  •  ,12,'  ,'  ,  12,  •  )  '  ) 
(I.J) 
E 
E 

A=3    1ATR  IX 
1=1 .NS 

,105)     I, J 

AX.'tNTEix    GA{  •  ,12,  «,  '  ,  12,  •  )'  ) 

(I.J) 

E 

HE    IDENT  ITY    MATRIX 

l.NS 

1,NS 

C.J)     GC    TO    300 

=0. 

,J)  =  0  . 

0. 

=1. 

,J)  =  1  . 

1. 

E 

E 


K.CT.NS)     GO    TO    13 


INPUT    OF    OEilREC    PULEo    LOCATION 


REAL    PART     OF     DESIRED    PCLE    LOCATION.*) 


IMAGINARY     PART    OF    DcSIaED    POLE    LOCATION.' I 


C 
C 


wRITE     (6,43) 

43  FORMAT {  •    ENTER 
READ, A 
WRITEI6  ,44) 

44  FORMAT*  •    ENTER 
READ,o 

BUILDING    ALPHA  (PHI) 

IF    (B     .EJ.OI     GO    TO     14 
Al=-2.*EX?{  a*T)  *COS(S*T) 
A2=EXP ( 2.*A*T ) 
EVALUATE    Pt-1»PHI 

CALL    VMULFF      ( PH ,PH ,NS , NS , NS , 20 ,20 , T M ,2 0, I ER ) 


601 

6C0 


DO    2     I  =  1,NS 
DJ    6    J=  1,NS 

TO(I  . J)  =  PH(  I, j)*Al 

T I  I  I , J)=AI  (  I, J)*A2 

TS(  I  ,  J)  =  TM(  I,  J  )-cTU 
CONT  INUE 

>,TS,NS,NS,NS,2  0,2  0,ATEHP,20, 


tli'l  ',  ji=Al  (  I,' J)*A2 

""II  I, J  J+TI  (  I  ,J) 


CCNf  1NUE 

CALL    VMULFF      ( ALPHA 

DO    600     1=1, NS 

DO    601     J=1,NS 

ALPHA(  I  ,J)  =  ATEMP(I  ,J) 

CONTINUE 


IfcR) 


ACK00C10 
AC  K.OoO  20 
ACK0CO3G 
AC\000<O 
ACKu0O50 
ACK.OJ060 
ACK.0G0  /o 
AC  KuuO  SO 
ACK0C090 
ACK00100 
ACK.0CU0 
ACK.CC120 
ACK0ol30 
ACKOOltO 
ACi<>001  50 
uCKO-l 60 
AC*00170 
ACKOOIBO 
AC*Oul<yO 
ACN00200 
AC*00210 
ACK0u220 
4C\0C230 
ACK0o240 
ACK00250 
ACKC3260 
AC  isG^2 70 
ACK.0J230 
AC  KG 02  90 
ACK0J3O0 
ACXGD310 
ACK00320 
ACn0j330 
ACK00340 
ACK00350 
AC  ,\0G360 
ACK003  70 
ACRC038G 
AC*00390 
ACK0O400 
ACK.0G4  1G 
ACKC0420 
ACK.DG430 
AC/vOD4^0 
ACKG0450 
AC  K0J46O 
a:\00470 
AC X0j4 80 
ACKO  J4  90 
ACK0C500 
ACK.00510 
ACK00520 
ACK0O530 
AC*0G540 
ACK00550 
ACK00560 
AC*DU570 
ACK0U580 
ACK00590 
ACK0U60D 
ACK00610 
ACR00620 
AC  K.0o630 
AC  K0U64O 
ACK00o50 
ACK00660 
AC  N0U6  70 
ACKOO60G 
AC  K00690 
ACK00700 
ACK00710 
ACKOU720 
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FILE:     ACKERMAN    wATFIV 


Ai 


14 


7 
3 


701 
700 


18 

9 

22 

10 


28 

11 


C 
C 

45 
SENTRY 


K=    K+2 

GO    TO     16 

Al=EXPl  A*T) 

DO    3     1  =  1,1 MS 

00    7    J=1,NS 

T0(  I  ,J)=Al (  I  ,  J) *Al 

TllltJl =PHl I , JJ-TO( I, J) 

CCNTINUt 

CALL    VMULFF  I  ALPHA,  TI,  NS, NS,N S  , 20 , 20 , AT EMP , 2 0,  IE?.  J 

DO    700     1  =  1, ,4S 

00    701     J=1,NS 

ALPHAi  I  ,JJ=ATEMr>(I  ,J  J 

CONTINUE 

K=K  +  1 

GO    TO    16 

00    9    1  =  1, NS 

J  =  l 

E(  I  ,  J»  =  GA(  I  ,  J  ) 

K=l 

IF    [K.GT.NS)     GC    TG    28 

00    10    1=1, NS 

CU  ,KJ  =  EI  i,l  ) 

CONTINUE 

K=K  +  1 

CALL    VMULFF     I PH , E, NS, NS , 1     ,  2  0  ,20  ,  E  ,  20,  IER  ) 

GO    TG    22 

DO    11    J=1,NS 

El     l.J)=AI I  NS,J  ) 

R=0. 

CALL    LGINF     (  C  ,2  0  ,.",S  ,  nS  ,R ,  D  ,  2  0  ,S  ,  «K  ,  I  Ef<  ) 


CALL    VMULFF 
CALL    VMULFF 


IE, 0,1     ,NS,NS,20,20,3F ,2G, IER) 

(HF  ,aL?HA,  1     ,  NS.:.S,2  0,20,UK,20,  IERJ 


PRINT,      •  CONTRUL    GAIN       VECTO* 

WRITE     (6,45)        I  OKI  1,J  J  ,J=1,NS  ) 

FORMAT  (  6F15  .4  ) 

RETURN 

END 


a:  \o 

ACkO 

ACKO 
ACkO 
ACnO 

ac.\o 

ACKO 

Au  K  G 
ACKO 
ACKO 
ACKO 
ACKO 
ACKO 
ACKO 
ACKO 
ACKO 
ACKO 
ACKO 
ACKO 
ACKO 
ACKO 
AC  KO 
aC/sO 
ACKO 
ACKO 
ACkO 
ACkO 
ACkO 
ACKO 
ACKO 
ACKO 
ACnO 
ACkO 


KO 
.KG 
,K0 
.KG 
,K0 
1K0 
,K0 


0  7  30 
\j  7  tO 
0  7  50 
u  7oO 
0770 
0730 
0790 
03  00 
Ca  10 
CS20 
0830 
o3  40 
Od50 
0360 
08  70 

08  8  0 
0390 
C900 

09  10 
0920 
j9  30 
0940 
0950 
G9t>0 
0970 
0980 
u990 
10C0 

10  10 
1020 
1030 
10^0 
1050 
1060 
1070 
10  30 
10  90 
1100 
lilO 
1120 
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APPENDIX  I 
COUPLED  STATE-FEEDBACK  DESIGN  INPUT  DATA  FOR  POPLAR  PROGRAM 

FILE:     FEEDBACK    OATA  Al 

1000  0 

200.0  0-0                     1.0 

1  0                         0 

0.4  0.4                   0.0  0 

0           0            0  5           3            3    LOCO            0 
17171702171702171  7 

1.5  /957586-01  3 .6 3693266-02  5.  58435  925" -03  -3 . 54  I  52 36E -2  -5 . 63 50634E -3 

1.4376905E-03  -o. 35oGl 9 5E- 3  3.12150  380  6+0  -2 .  743h69oE +0  3 .496u341 6 -01 

0  0  0  0  0 

0  0 

7.3702121E-03  1  .0001 3475-00  1  .  2  L45  7  1  35 -02  -1.3  110  3295-4  -2 . 3707903E-5 

6.48032775-06  -4 . 3 14336 7E-5  2. 35 544 1 5c -02  -2.16  5  79256-2  9.078*7435-03 

0  0  0  0  3 

0  0 

-3.5755072E-2  -8. 968 3556E-4  9 . 3927 5 54E -0 1  8  .  7 9 282 40E- 04  1 .3 9 377 5 5c-04 

-3.  1463045E-5  2.094o99bc-J4  -  1. 1 434246 6- 1  1 .05 134 94c- 01  -4 .4066 55 1 c-2 

0  0  0  0  0 

0  0 

6.3641914E-04  4 . 3072347E-03  7 . 37304 52E -04  9 . 1 1  345 92 fc-0  1  1 . 1 7240 94c-02 

1.03891730E-7  -  1. 02 76906E-O  -  2.  5  7  72  7  _J3  E-  1  -4  .  3  7  j>06  7o  5  -  1  3.  o2396946-04 

0  0  0  0  0 

0  0 

-1.2632605E-4  -6. 3  10 21 33 E-4  -  1. J508452 5-4  -5.  5242591 E-l  9.  94o 1 43 15-01 

-1.1935039E-3  1  .34677576-0 7  1 . 1052479E -01  1 . 0373553E- 01  -5 .65300o2 E-5 

0  0  0  0  0 

0  0 

5.0C58503E-05  3.4906 702E-04  5 . 3521 3 3 5E -05  -  3 . 42 3525 76 -4  -5 . 4232 1 94P-5 

9.9653837E-01  1.2  473  8  785-02  -4.3337960 E-2 '  -4. 3 103 149E -2  2 .2 3902230E-5 

0  0  0  0  0 

0  0 

1. 1198072E-02  8  .45 8971 5E-02  1  .298 26 64E -02  -  3 . 2 578 336E-2  -1 . 3  I 94430E-2 

-5.53283155-1  9  .9465984E-01  -4. 5 3 1 3 1 72 E-0  -  7 . 69 330 73 E -0  6.37744  246-03 

0  0  0  0  0 

0  0 

-1.8197865E-3  -  1. 6475936E- 2  -  2. 5342 871 E-3  1 . 59 144 37E- 02  2.592  74  895-03 

-7.69372665-7  o  .  1 1292605-06  1  .  9q 582 56b -02  3.25513265-01  -1.61176915-3 

0  0  0  0  0 

0  0 

-2. 4201  113E-3  -2.62507035-2  -4.34o3 1 1 1 5- 3  2.5  073  7845-02  4. 1 6 749 755-0 i 

-2.20856125-6  1  .4o5903*6-05  -  1.  0  1  J4  5  58  5-  I  8.5  2133  305-01  -3  .  0o06-4  3  75- 3 

0  0  0  0  0 

0  0 

8.51203615-03  7 . 7  170 2C6E-02  1  .  18702 99E -02  -7 . 45350666 -2  - 1 .2 1 44650E -2 

3.46475185-03  -1.25026405-2  4  .65o J4625 -Ob  -3 . 84o95 975 -0  I . Oo 756 1 05-00 

0  0  0  0  0 

0  0 

0  0  0  0  0 

n  "l  "i  o  ^ 

9.97738  786-01  -2.  7CI 3690E-3  1.02458  966-03  2.2  68  5  7365-01  -2. 2  392  473E-1 

3.2285155E-C0  3  .u730433E-02 

0  0  0  0  0 

0  0  0  0  0 

1.249  2612E-02  9.9999  1215-01  3.36863  395-06  7  .  92  133  8  IE- 04  -7.  325  1  0  1  9E-4 

1. 61206056-02  3.2l4o5656-04 

0  0  0  0  0 

0  0  0  0  0 

-1.32996045-3  -2.093  57296-1  9.04337  2  76-01  -3.68319656-5  3.64511286-05 

-9.65973186-4  -2.6  9033  705-5 

0  0  0  0  0 

0  0  0  0  0 

-9.07977235-3  -  I. 102408 7E-2  3  .9997292E -03  9. 2 763032E- 0  1  1 . 1 630732E-02 

-1.70079726-1  -2. 0028512E-3 

0  0  0  0  0 

0  0  0  0  0 

-1. 7340942E-4  -2.0717899E-4  7 . 85645 38E -05  1 . 7 374049E-02  9. 1 0594 3 1 E-Ol 

2. 45055106-01  2.72575106-03 

0  0  0  0  0 

0  0  0  0  0 

-1.90810306-3  -2. 30o 18305-3  8.4o098055 -04  1.5  2  307156-01  -1.50  374186-1 

7.46927475-01  -1.142948  56-3 

0  0  0  0  0 

0  0  0  0  0 
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FILE:     FEEDBACK    JATA 


Al 


-b. 28757186-2 

- 

7. 56C2497E-I 

J..817ol9l6 

-0  2 

5.oi224356- 

-00 

-5.531^0555-0 

4.33236486*01 

5 

.1J4J68  56- 

■02 

-7.56093296-3 

0 

-1.64491866-2 

0 

7.98216236-02 

0 

-9.  9881  6086 -<t 

0 

1.8285612E-04 

0 

-7.2448  7266-5 

0 

-1.75755156-2 

0 

3.«t3l55566-C3 

0 

5.43075556-C3 

0 

-l.o0730196-2 

0 

0 

2 

.70136906- 

-33 

0 

d 

.792o971c- 

■0  6 

0 

2 

.09357296- 

-Gl 

0 

1 

.10240ti75- 

-0  2 

0 

2 

.0717399E- 

-0  4 

0 

2 

.3081  8306- 

-0  3 

0 

7 

.5oQ24976- 

-0  2 

1.00C00006-C0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1 

.ocoooooe- 

-JO 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1  .OOOOOOOE 

-00 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1 .00000006- 

-00 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1.00003 006-00 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

I.  0  COCO  006-00 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

I 

.GCCU0006- 

-30 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1 .OOOOOOOE 

-00 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1  .0000000c- 

-00 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1.00003  006-00 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1.0000000E-00 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1 

.OOOOOQOE- 

-30 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1 .00000006 

-00 

0 

0 

0 

0 

0 

0 

0 

0 

0 
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FILE:     FEEDBACK.    OA  T  A 


Al 


1. OOOOO OCE 


-0.0 

0.0 

0.0 

0.0 

O.J 

CO 

CO 

-1.4 

0    0    0 

0    0    0 

-1 

0.13 
0.13 
0.89 
0.b9 
0.95 
0.95 
0.99 
1.00 
0.59 
0.99 
0.09 
0.33 
0.83 
0.96 
0.9o 
0.39 
0.93 


188000E 
OOOOOOE 
000000c 
OOOOOOE 
OOOOOOfc 
OOOOOOt 
COOOOOt 
S90000E 
0    0    0 
0    0    0 
0000 
199300 
199300 
698700 
698700 
362500 
362500 
8202  00 
000000 
560800 
360800 
92  80  5  0 
521600 
321600 
9C8700 
908700 
379700 
707500 


0 
0 
0 
0 
0 
0 
0 
0 
0 
0 

-co 

0 

0 

0 

0 

-00 

♦  00 
+  00 
+  C0 

♦  00 
-GO 

♦  CO 
-00 

0    0    0 

coo 

100 
0.03 

-COi 
G.09 

-C09 
C.02 

-C  02 
0.00 
COO 
0.03 

-COB 
COO 
C  jO 

-C.30 
COi 

-C.03 
COO 

coo 


l.OCOO 
-0.0279 
O.OCOl 
0  .  0  C  00 
0.000 
O.OCCJ 

o.occo 

0  .  OGG-» 
0.0C23 
0    0    0 
0   0    c 
00 

193770 
153770 
067030 
o67080 
"J  9  85  00 
S9S500 
000000 
JCCOOO 
307130 
307130 
OOOOuO 
573600 
v  7  36  0  0 
235330 
235830 

oocooo 
occooo 


0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

000E-J3 

OOOc-CO 

OOOE+00 

OJOE+OO 

OOOOE+0 

OOOE+OQ 

OGGc+Ou 

OOOiz-00 

OOOE-00 

0    0    0    0 

0    0    0    0 


0 
0 

0 
0 
0 

0 
0 
0 

0 
0 
0 

0.0042300E-00 
1.7732000E+-00 
O.OOOOOObt  +-00 

O.COOGOOOE+OO 

C.0003GOOE+0 

0.G00G3Q0E+GG 


1.03GJOO0E-u0 

0 
0 
0 

0 
0 
0 

0 
0 
0 

0.  02  900  00 1:-*- 00 

■2.016200Gc  +  00 

G.GGoOGOOc+OO 

0.00000006+0 
0. OOOOOO^E+O 
0. 01 78000E-0 


0 

0 

0 

0 

1.0000300E-00 

0 
0 
0 

0 
0 
0 

0.00413C0E+00 

o.uii ioooe+oo 

0.00030QGt*00 

0. 30033 OOE+OO 

0 . OGuOOOOE+0 

-0. Jl 79000E-0 
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APPENDIX  J 
COUPLED  ESTIMATOR  DESIGN  INPUT  DATA  FOR  POPLAR  PROGRAM 


FILE:     ESTIMATQ    J«TA 


Al 


10CO  0 

100.0 

1 

0.4 

0  0  0 

17171702171702  1717 

1.5795758E-01 

1.43769U5E-03 

0 

0 

7.3702121E-03 

6.4803277E-C6 

0 

0 

-3.5755072E-2 

-3. 1463045E-5 

0 

0 

6.3641914E-G4 

1.03891780E-7 

0 

0 

-1.2632  605E-4 

-I. 193503 9E-8 

0 

0 

5.0C58503E-05 

9.9653887E-01 

0 

0 

1.1 1980  72E-C2 

-5.5328815E-1 

0 

0 

-1. 8197865E-3 

-7.o93726dE-7 

0 

0 

-2.4201 118E-3 

-2.2035612E-6 

0 

0 

8.5120361E-03 

3.464751ot-03 

0 

0 

0 

0 

9.97738  78E-C1 

3.2285155E-C0 

0 

0 

1.2492612E-02 

1.6120605E-02 

0 

0 

-1.3299604E-3 

-9.6597318E-4 

0 

0 

-9.C797723E-3 

-1.  7007972E-1 

0 

0 

-1.7340  942E-4 

2.4505510E-Q1 

0 

0 

-1.9081030E-3 

7.4692747E-01 

0 

0 


0.0 
G 

0.4 
5 


1.0 
0 

0.0 
3 


0 
1000 


3.63c932bc-J2 

-o.  3550195E-3 

0 

0 

1  .0001 347E-00 
-4. J148367E-5 

0 

0 

-e.  9b3355bE-4 

2  .0-^46  99bt-04 

0 

0 

4.8072347E-03 

-1.  02  7690oE-b 

0 

0 

-8.S  1021d3E-4 

1.J487757E-J7 

0 

0 

3.49G6702E-04 

1.2473878E-02 

0 

0 

8.4589715E-02 

9.9465  9b<»E-0  1 

0 

0 

-1.6475936E-2 

6.1 1292o0E-06 

0 

0 

-2.6250703E- 2 

1.ho59039E-05 

0 

0 

7.7170206E-02 

-1.250<.b40E-2 

0 

0 

0 

0 

-2. 7013690E-3 

3  .6730-O3E-J2 

0 

0 

9.9999121E-01 

3.214o565E-u4 

0 

0 

-2.0935729E-1 

-2.69J8370E-5 

0 

0 

-1.  1024087E-2 

-2.0028512E-3 

0 

0 

-2.  071  7899E-4 

2.7257510E-03 

0 

0 

-2.3081830E-3 

-1.  1429485E-3 

0 

0 


5.5J^d592E-03 

3.1215G380E+0 

0 

1.2145713E-02 
2. J5 544 15E -02 

u 

9.3927554E-01 
-1.143424oE-l 

0 

7.37e0452E-04 

-2. 5772733 E- I 

0 

-1.350d452E-4 

1 .1052479E-01 

0 

5.3521335E-05 

-4.3d37960E-2 

0 

1.2932664E-02 

-h. 53131  7<.E-G 

0 

-2.5342371E-3 
1  .  9b5325bE-02 


-4.04031  11  E-3 

-1.0134554E-1 

0 

1  .  1870299E-02 

4  ,65b:><t62E-00 

0 

0 

0 

1.0245896E-03 

0 

0 

3.36  8b339E-06 

0 

0 

9.0483727E-01 

0 

0 


-3.541523oE-2 
-2. 7434696E+0 

0 

-i. 3110829E-4 

-2.1657925E-2 

0 

8.7523240E-04 
1.051 J494E-01 

0 

9.9184592E-01 

-4. 3730676E-1 

0 

-5.  5242591E-1 

1.0t37355dt-01 

0 

-3.4235257E-4 
-t. 3103149E-2 

0 

-3.2578335E-2 

-7.6933073E-0 

0 

1. 5914437E-02 
d.2551326E-Cl 


2.507373-*E-02 

d.5213330c-01 

0 

-7.'»535066E-2 

-3.o409597E-0 

0 

0 

0 

2.2685736E-01 

0 

0 

7.9213831E-04 

0 

0 

-3.68319b5E-5 

0 

0 


3.9997292E-03  9 . 2 7b30 32E-0 1 


0 
0 


7.8564538E-05  1  .  7 374049E -02 


0 

0 

8.  46  098  0  5c -04 

0 
0 


0 

0 

1.5280715E-01 

0 
0 


-5.b350634E-3 
d.-»9b0341c-ul 

0 

-2.3707903E-5 

9.0789743E-03 

0 

1.3937755E-04 

-t.40bb  551E-2 

0 

1.1724094E-02 

3.o239694t-C4 

0 

9.94iil481E-01 
-5.6330062E-5 

0 

-5.4232194E-5 
2.23902230E-5 

0 

-1.3194430E-2 

6.3774424E-03 

0 

2.5927439E-03 

-1.611 769lt-3 

0 

4. 1674975E-03 

-3.0605437E-3 

0 

-1.2144650E-2 

1.0075510E-00 
0 

0 

0 

-2.2J92478E-1 

0 

0 

-7.8251 019E-4 

0 

0 

3.6451128E-05 

0 

0 

1.1630732E-02 

0 

0 

9.1059431E-01 

0 

0 

-1.5037418E-1 

0 

0 
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FILE:     ESTIMATJ    DATm 


Al 


-6.2875713E-2 
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